2105.13549v3 [astro-ph.CO] 21 Mar 2022 


e 
. 


arXiv 


DES-2020-0617 
FERMILAB-PUB-21-221-AE 


Dark Energy Survey Year 3 Results: 
Cosmological Constraints from Galaxy Clustering and Weak Lensing 


T. M. C. Abbott,! M. Aguena,” A. Alarcon,’ S. Allam,* O. Alves, 5? A. Amon,’ F. Andrade-Oliveira,®? J. Annis,* S. Avila, 
D. Bacon,’ E. Baxter,!° K. Bechtol,!! M. R. Becker,? G. M. Bernstein,!* S. Bhargava,!? S. Birrer,!^ J. Blazek,'^:!ó 
A. Brandao-Souza,!^? S. L. Bridle,!® D. Brooks? E. Buckley-Geer,?^ ^ D. L. Burke,^?! H. Camacho,®? A. Campos,” 
A. Carnero Rosell,?:??^ M, Carrasco Kind,?°»° J. Carretero,?" F. J. Castander,?2*?? R. Cawthon,!! C. Chang,”” 30 A. Chen? 
R. Chen,?! A. Choi,? C. Conselice,!*:33 J. Cordero,!® M. Costanzi, 43536 M. Crocce,”*:29 L. N. da Costa? ?" M. E. da 
Silva Pereira? C. Davis,’ T. M. Davis,?? J. De Vicente,? J. DeRose,® S. Desai,*! E. Di Valentino,!'® H. T. Diehl,* 
J. P. Dietrich, S. Dodelson,?^ 9 P. Doel,? C. Doux," A. Drlica-Wagner,?? ^? K, Eckert,!? T. F. Eifler,^^ ^ F. Elsner,!? 
J. Elvin-Poole,??49 S. Everett,” A. E. Evrard, ^9? X. Fang, ^ A. Farahi,^ ^? E. Fernandez, I. Ferrero,?? A. Ferté,“ 
P. Fosalba,??:?? O. Friedrich?! J. Frieman,**° J. Garcia-Bellido,® M. Gatti,'? E. Gaztanaga, ^5 7? D. W. Gerdes,^?? 

T. Giannantonio,°” >! G. Giannini," D. Gruen;?? ^?! R. A. Gruendl,??6 J. Gschwend,? 3” G. Gutierrez,^ I. Harrison,?^ !8 
W.G. Hartley,?? K. Herner,* S. R. Hinton? D. L. Hollowood,*” K. Honscheid,?- 46 B. Hoyle,*” E. M. Huff,? D. Huterer,> 
B. Jain, D. J. James;?? M. Jarvis, N. Jeffrey, !” 57 T. Jeltema,” A. Kovacs,?^?^ E. Krause,“ R. Kron,^ ? K. Kuehn,’ >? 

N. Kuropatkin,‘ O. Lahav,!? P.-F. Leget,’ P. Lemos,!™! A. R. Liddle, 61-62 C, Lidman,®?° M. Lima,? H. Lin,‘ 

N. MacCrann,® M. A. G. Maia,>3” J. L. Marshall, P. Martini? 95:6? J. McCullough,’ P. Melchior,’° J. Mena-Fernández,?? 
F. Menanteau,”>»*° R. Miquel,’ ?7 J. J. Mohr,?:? R. Morgan,!! J. Muir,’ J. Myles,?? ^?! S. Nadathur,!? A. Navarro-Alsina, 
R. C. Nichol,’ R. L. C. Ogando,^?7 Y. Omori,®3%7 A. Palmese,^?? S. Pandey,!? Y. Park,? F. Paz-Chinchón,?^ *? 

D. Petravick,? A. Pieres,??" A. A. Plazas Malagón,” A. Porredon,?? 46 J. Prat??? M. Raveri,!? M. Rodriguez-Monroy,”? 
R. P. Rollins,! A. K. Romer,!? A. Roodman,”! R. Rosenfeld," ^? A. J. Ross,*” E. S. Rykoff,^?! S. Samuroff,?? C. Sánchez," 
E. Sanchez,” J. Sanchez,* D. Sanchez Cid,’ V. Scarpine,* M. Schubnell,” D. Scolnic,3! L. F. Secco,/?:99 S, Serrano,?®: 79 
I. Sevilla-Noarbe,? E. Sheldon,” T. Shin, M. Smith," M. Soares-Santos,? E. Suchyta,"" M. E. C. Swanson,” M. Tabbutt,!! 
G. Tarle? D. Thomas,” C. To? 7?! A, Troja," ^? M. A. Troxel,?! D. L. Tucker,’ I. Tutusaus,?*?? T. N. Varga, ^78 
A. R. Walker,! N. Weaverdyck,? R. Wechsler,??: ^?! J. Weller, 7? B. Yanny,* B. Yin,? Y. Zhang,* and J. Zuntz?? 


(DES Collaboration) 


! Cerro Tololo Inter-American Observatory, NSF's National Optical-Infrared Astronomy Research Laboratory, Casilla 603, La Serena, Chile 
?Laboratório Interinstitucional de e-Astronomia - LIneA, 
Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil 
Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA 
4Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA 
5 Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 
Instituto de Física Teórica, Universidade Estadual Paulista, Sáo Paulo, Brazil 
Kavli Institute for Particle Astrophysics & Cosmology, 
P. O. Box 2450, Stanford University, Stanford, CA 94305, USA 
8 Instituto de Fisica Teorica UAM/CSIC, Universidad Autonoma de Madrid, 28049 Madrid, Spain 
? Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth, POI 3FX, UK 
V Institute for Astronomy, University of Hawai’i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 
" Physics Department, 2320 Chamberlin Hall, University of Wisconsin-Madison, 1150 University Avenue Madison, WI 53706-1390 
12 Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA 
P Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BNI 9QH, UK 
V Graduate School of Education, Stanford University, 160, 450 Serra Mall, Stanford, CA 94305, USA 
P Department of Physics, Northeastern University, Boston, MA 02115, USA 
16 Laboratory of Astrophysics, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland 
V Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, 13083-859, Campinas, SP, Brazil 
'8Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, 
University of Manchester, Oxford Road, Manchester, M13 9PL, UK 
P Department of Physics & Astronomy, University College London, Gower Street, London, WCIE 6BT, UK 
? Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA 
?! SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA 
? Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA 
°3 Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain 
" Universidad de La Laguna, Dpto. AstrofAsica, E-38206 La Laguna, Tenerife, Spain 
P Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA 
?6 Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 W. Green Street, Urbana, IL 61801, USA 
?' Institut de Física d'Altes Energies (IFAE), The Barcelona Institute of 
Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain 
?5 Institut d'Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain 


” Institute of Space Sciences (ICE, CSIC), Campus UAB, 
Carrer de Can Magrans, s/n, 08193 Barcelona, Spain 
9 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA 
?! Department of Physics, Duke University Durham, NC 27708, USA 
P Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA 
P University of Nottingham, School of Physics and Astronomy, Nottingham NG7 2RD, UK 
9 Astronomy Unit, Department of Physics, University of Trieste, via Tiepolo 11, I-34131 Trieste, Italy 
55 INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, I-34143 Trieste, Italy 
P6 Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy 
37 Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil 
38 School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 
? Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain 
Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 
^! Department of Physics, IIT Hyderabad, Kandi, Telangana 502285, India 
? Faculty of Physics, Ludwig-Maximilians-Universitdt, Scheinerstr. 1, 81679 Munich, Germany 
9 NSF AI Planning Institute for Physics of the Future, 
Carnegie Mellon University, Pittsburgh, PA 15213, USA 
^ Department of Astronomy/Steward Observatory, University of Arizona, 
933 North Cherry Avenue, Tucson, AZ 85721-0065, USA 
“Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA 
“© Department of Physics, The Ohio State University, Columbus, OH 43210, USA 
“7 Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064, USA 
5 Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA 
? Departments of Statistics and Data Science, University of Texas at Austin, Austin, TX 78757, USA 
P Institute of Theoretical Astrophysics, University of Oslo. P.O. Box 1029 Blindern, NO-0315 Oslo, Norway 
?! Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 OHA, UK 
? Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, UK 
3 Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA 
? Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK 
? Department of Astronomy, University of Geneva, ch. d'Écogia 16, CH-1290 Versoix, Switzerland 
"Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 
>! Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, 
Université PSL, CNRS, Sorbonne Université, Université de Paris, Paris, France 
58 Australian Astronomical Optics, Macquarie University, North Ryde, NSW 2113, Australia 
? Lowell Observatory, 1400 Mars Hill Rd, Flagstaff, AZ 86001, USA 
8 Institute for Astronomy, University of Edinburgh, Edinburgh EH9 3HJ, UK 
9! Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, 1769-016 Lisboa, Portugal 
8 Perimeter Institute for Theoretical Physics, 31 Caroline St. North, Waterloo, ON N2L 2Y5, Canada 
9 Centre for Gravitational Astrophysics, College of Science, 
The Australian National University, ACT 2601, Australia 
8 The Research School of Astronomy and Astrophysics, Australian National University, ACT 2601, Australia 
$5 Departamento de Física Matemática, Instituto de Física, 
Universidade de Sáo Paulo, CP 66318, Sáo Paulo, SP, 05314-970, Brazil 
$6 Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 OWA, UK 
© George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, 
and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77845, USA 
$88 Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA 
© Radcliffe Institute for Advanced Study, Harvard University, Cambridge, MA 02138 
"Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA 
7 Institució Catalana de Recerca i Estudis Avançats, E-08010 Barcelona, Spain 
® Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany 
P Kavli Institute for the Physics and Mathematics of the Universe (WPI), 
UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 
"ICTP South American Institute for Fundamental Research 
Instituto de Física Teórica, Universidade Estadual Paulista, São Paulo, Brazil 
P Brookhaven National Laboratory, Bldg 510, Upton, NY 11973, USA 
School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK 
™ Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 
18 Universitdts-Sternwarte, Fakultät für Physik, Ludwig-Maximilians Universitit München, Scheinerstr. 1, 81679 München, Germany 
(Dated: March 22, 2022) 


We present the first cosmology results from large-scale structure using the full 5000 deg? of imaging data 
from the Dark Energy Survey (DES) Data Release 1. We perform an analysis of large-scale structure combining 
three two-point correlation functions (3x2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy 
clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions, galaxy—galaxy lens- 
ing. To achieve the cosmological precision enabled by these measurements has required updates to nearly 
every part of the analysis from DES Year 1, including the use of two independent galaxy clustering samples, 
modeling advances, and several novel improvements in the calibration of gravitational shear and photometric 
redshift inference. The analysis was performed under strict conditions to mitigate confirmation or observer bias; 
we describe specific changes made to the lens galaxy sample following unblinding of the results and tests of 
the robustness of our results to this decision. We model the data within the flat ACDM and wCDM cosmo- 
logical models, marginalizing over 25 nuisance parameters. We find consistent cosmological results between 


the three two-point correlation functions; their combination yields clustering amplitude Sg = pru n 
and matter density Qm = 0.339*0022 in ACDM, mean with 68% confidence limits; Ss = 0.7750 025, 
Qn = (52733 and dark energy equation-of-state parameter w = —0.98* 0:32 in wCDM. These con- 


straints correspond to an improvement in signal-to-noise of the DES Year 3 3x2pt data relative to DES Year 
1 by a factor of 2.1, about 20% more than expected from the increase in observing area alone. This combina- 
tion of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave 
background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed p — 0.13 to 0.48. 
We find better agreement between DES 3x2pt and Planck than in DES Y1, despite the significantly improved 
precision of both. When combining DES 3 x2pt data with available baryon acoustic oscillation, redshift-space 
distortion, and type Ia supernovae data, we find p — 0.34. Combining all of these data sets with Planck CMB 


lensing yields joint parameter constraints of Sg = 0.812908, Qn = 0.306+9:005, h= 0.680+9:003, and 


Sm, < 0.13 eV (95% CL) in ACDM; Sg = 0.81 


w = —1.031*0022 in wCDM. 


I. INTRODUCTION 


The discovery of the accelerated expansion of the universe 
led to a new standard model of cosmology, which 
is dominated by a spatially smooth component with nega- 
tive pressure called dark energy. Over the intervening two 
decades, the evidence for the presence of dark energy has be- 
come much stronger thanks to data from an impressive variety 
of cosmological probes. Modern cosmological measurements 
using type Ia supernovae [BHIT], cosmic microwave back- 
ground (CMB) fluctuations [121/14], galaxy clustering 
, and weak gravitational lensing are in agreement 
with a spatially flat universe with about 30% matter (visible 
and dark) and 70% dark energy. 


However, the physical nature of the dark energy that causes 
accelerated expansion remains unknown. The simplest and 
best-known phenomenological model for dark energy is the 
energy density of the vacuum, incorporated in the field equa- 
tions of General Relativity by the cosmological-constant term 
A [29]. The resulting A Cold Dark Matter (ACDM) model 
serves as a benchmark for tests with current and future data. 
Beyond ACDM, there exists a rich set of other potential 
models to explain cosmic acceleration, including evolving 
scalar fields, modifications to general relativity, and other 
physically-motivated possibilities. This has spawned an ac- 
tive research area focused on describing and modeling dark 
energy and its effects on the expansion rate and the growth of 


density fluctuations [30][31]. 


The quest to understand dark energy has spawned a world- 
wide effort to better measure the growth and evolution of 
cosmic structure in the universe. The current generation of 
observations is spearheaded by the so-called Stage-III dark 
energy experiments, which include the Dark Energy Sur- 
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vey (DES! [32134], the Hyper Suprime-Cam Subaru Strate- 
gic Program (HSC [36], the Kilo-Degree Survey 
(KiDS]"|[28] [37], and the Extended Baryon Oscillation Spec- 
troscopic Survey (eBOSS}'] (38). These surveys have demon- 
strated the feasibility of ambitious large-scale structure anal- 
yses, and featured extensive tests of theory, development of 
state-of-the-art systematics calibration, and new rigor in pro- 
tecting analyses against observer bias before the results are 
revealed. These surveys have, thus far, provided constraints 
consistent with the ACDM model, and contributed to tighten- 
ing the constraints on several of the key cosmological param- 
eters related to dark matter and dark energy. 

Large-scale structure (LSS) in the universe provides a 
powerful set of tools to probe dark energy. The statistics 
and temporal growth of cosmic structure complement the 
largely geometrical sensitivity to dark energy of type Ia su- 
pernovae and the CMB. For nearly half a century, measure- 
ments of the galaxy two-point correlation function, a statis- 
tic describing the spatial clustering of galaxies, have provided 
pioneering cosmological constraints and early evidence for 
the ACDM model [391/52], as well as recent, high- 
precision constraints on the cosmological parameters 
. Another aspect of LSS that is sensitive to both dark mat- 
ter and dark energy is cosmic shear, slight distortions of the 
shapes of distant background galaxies due to weak gravita- 
tional lensing of light passing through the structures between 
these sources and us. While the interpretation of galaxy clus- 
tering is complicated by galaxy bias [67], cosmic shear 
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measurements are more directly related to the distribution of 
mass. First detections of cosmic shear have been 
followed by an impressive maturing of this probe, with in- 
creasingly more competitive constraints on cosmological pa- 
rameters [23125] (72H78]. Finally, galaxy—galaxy lensing, the 
cross-correlation of lens galaxy positions and source galaxy 
shapes, provides a link between galaxy clustering and cosmic 
shear. Galaxy-galaxy lensing measurements have also ma- 
tured to the point where their combination with galaxy clus- 
tering breaks degeneracies between the cosmological param- 
eters and bias, thereby helping to constrain dark energy 
[92]. The combination of galaxy clustering, cosmic shear, and 
galaxy-galaxy lensing measurements powerfully constrains 
structure formation in the late universe, while strongly self- 
calibrating many astrophysical or systematic parameters in the 
model. 


The stakes have become higher with recent evidence for 
possible tensions between parameters as measured by dif- 
ferent cosmological probes. These tensions may indicate 
new physics beyond ACDM — or else could be due to 
unaccounted-for systematics or an underestimation of uncer- 
tainty in some probes. Potentially most significant among 
these is the “Hubble tension," indicated by a ~4-60 dis- 
crepancy between measurements of the Hubble constant in- 
ferred from the primary CMB anisotropies and higher 
values measured from a local distance anchor, such as, most 
prominently the astronomical distance ladder (e.g., [4p 
or masers [95], though some measurements also indicate a 
lower value in better agreement with the CMB. The 
Hubble tension may indicate new physics, and it is crucial to 
improve measurements, revisit assumptions and systematics 
[e.g., and invest in novel, independent methods and 
probes [99] [100]. 


Additionally, several experiments that are sensitive to the 
growth of structure have historically preferred, on average, 
lower values of the parameter Sg = og(Qm/0.3)°° rela- 
tive to that predicted by the CMB anisotropy, where the am- 
plitude of mass fluctuations og is scaled by the square root 
of matter density Qm. This parameter is predicted to be 
higher with the CMB than is measured in lensing (e.g., 
(271[33][101][102]). The difference has been claimed by other 
experiments to be as large as 2-30. Other probes of the late 
universe, in particular spectroscopic galaxy clustering [103], 
redshift-space distortions (RSD) [21], and the abundance of 
galaxy clusters [104] [105], also tend to favor a lower Sg than 
that measured by the CMB on average (assuming the ACDM 
model). 

Previously, the DES Collaboration analyzed data from its 
first year of observations, which covered 1514 deg?, and 
constrained cosmological parameters using galaxy clustering 
and gravitational lensing in ACDM and wCDM [33], con- 
strained beyond-wCDM models [106], and carried out nu- 
merous other tests of the standard cosmological framework 
[107/114]. Along with the aforementioned 
KiDS and HSC observations and analyses, the DES Y1 anal- 
ysis emphasized redundancy using two shape measurement 
methods that are independently calibrated, several photomet- 
ric redshift estimation and validation techniques, and two in- 
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dependent codes for predicting the measurements and per- 
forming a likelihood analysis. 


This paper presents key cosmological constraints from the 
first three years of observations (henceforth Y3) of DES. The 
DES Y3 data set analyzed here uses images covering nearly 
5000 sq. deg., or more than three times the area of Y1. It also 
dramatically increases the number of source and lens galaxies, 
and introduces new techniques for the analysis and treatment 
of statistical and systematic errors. As in Y1, we rely on a 
key cosmological probe of photometric LSS surveys, the so- 
called ‘3x2pt’ analysis, consisting of three two-point corre- 
lation functions: (i) w(0), the angular correlation function of 
the lens galaxies; (ii) y;(@), the correlation of the tangential 
shear of sources with lens galaxy positions; and (iii) €+ (0), 
the correlation functions of different components of the ellip- 
ticities of the source galaxies. We use these measurements 
only on large angular scales, for which we have verified that 
a relatively simple model describes the data, although even 
with this restriction we must introduce 25 free parameters to 
capture astrophysical and measurement-related systematic un- 
certainties. The paper is built upon and uses tools and results 
from 29 accompanying papers that are summarized 
in App. [A] We summarize in App. [B]|the major updates to the 
analysis that are different from the DES Y1 3x2pt analysis. 


The cosmological quantity that is best constrained by the 
3x2pt analysis is the overall amplitude of matter clustering in 
the low redshift universe, parameterized by Sg. The precise 
measurement of Sg in this paper allows a powerful test for 
consistency between the growth of structure and the expan- 
sion history in the broad class of cosmic acceleration models 
based on General Relativity (GR) and dark energy. Imple- 
menting this test requires a CMB anchor for the matter clus- 
tering amplitude at high redshift, and the test becomes sharper 
and more general when supernova and baryon acoustic oscil- 
lation (BAO) data are used to constrain the expansion history. 
DES probes matter clustering out to z œ% 1, so it also con- 
strains dark energy models on its own through the history of 
structure growth over this redshift range. The degeneracy be- 
tween Qm and gg in Sg is broken partly by this redshift evolu- 
tion and partly by the shape of the correlation functions, and 
it can be broken more strongly using external data that are 
sensitive to Qm. Lensing measurements depend on the ex- 
pansion history through the distance-redshift relation. This 
dependence affects our analysis, but the geometric constraints 
from DES weak lensing are not as strong as those from current 
supernova and BAO data. 


In subsequent sections of the paper, we focus first on the 
DES data sets in Sec. [II]and measurements of the three two- 
point correlation functions in Sec. We describe the mod- 
eling and analysis in Sec. then turn to the primary re- 
sults, tests, and parameter constraints from combining these 
measurements with additional measurements from DES, the 
CMB, and other external supernova, BAO, and RSD data in 
Sec.[V] We conclude in Sec. [VIII] 


Il. DARK ENERGY SURVEY DATA 


The Dark Energy Survey (DES) was a six-year observing 
program using the Dark Energy Camera (DECam [32]) at 
the Blanco 4m telescope at the Cerro Tololo Inter-American 
Observatory (CTIO) in Chile. The survey covered 5000 
deg? in grizY bandpasses with approximately 10 overlap- 
ping dithered exposures in each filter (90 sec in griz, 45 sec 
in Y) covering the survey footprint. In this paper, we uti- 
lize data taken during the first three years of DES operations 
(DES Y3), which made up DES Data Release 1 (DR1 [144]). 
This analysis uses imaging data covering the full 5000 deg? 
survey footprint for the first time, at approximately half the 
full-survey integrated exposure time. Preparing the imaging 
data for cosmological analysis is an exacting, multi-year pro- 
cess, and analysis of the final six-year data set is now in its 
early stages. The data is processed, calibrated, and coad- 
ded to produce a photometric data set of 399 million objects 
that is further refined to a ‘Gold’ sample for cosmological use 
115] [145] [146]. The Gold sample includes selection require- 
ments (cuts) on minimal image depth and quality, additional 
calibration and deblending, and quality flags to identify prob- 
lematic photometry and regions of the sky with substantial 
photometric degradation (e.g., around bright stars). The Gold 
galaxy sample extends to a signal-to-noise > 10 (extended) 
limiting magnitude of 23 in i-band. The final Gold sample 
used in this work after all cuts contains 319 million objects. 


In addition to the wide-field Gold sample, we rely on data 
from the DES deep fields covering a subset of the 27 
deg? DES transient search regions and the separate COSMOS 
field [147]. These images are taken in ugrizY bandpasses 
with DECam, and also have overlapping VIDEO or Ul- 
traVISTA imaging for near-IR photometry in Y JH K 
over most of the area. Coadd images are constructed from the 
best images (i.e., with smallest point-spread function (PSF) 
full-width half-maximum (FWHM)) with a goal of attaining a 
depth approximately 10x the typical wide-field coadd image 
depth. From these coadd images, we produce a deep catalog 
of 2.8 million objects that has 10c limiting magnitude of 25 
in i-band, and photometric variance 0.1 x the typical wide- 
field variance. This catalog helps to validate and calibrate 
our wide-field data in several ways. First, it is used to create 
an input model space for representative objects to draw onto 
wide-field-quality images in the Balrog or weak lens- 
ing image simulations [120]. Balrog is used to test the survey 
selection function, which describes the probability that an ob- 
ject type drawn from a complete galaxy population will be 
detected in our wide-field survey, and the weak lensing sim- 
ulations are used to test our shear calibration. The deep cat- 
alog also serves as a stepping-stone in our redshift inference 
methodology [123]. It allows us to map available spectro- 
scopic or many-band deep photometric observations into the 
ugrizJ HK bandpass space of our deep catalog, for which 
we have 1.68 million sources with matched near-infrared pho- 
tometry covering an area of 5.88 deg?. This is then mapped 
through Balrog onto wide-field galaxy information. 


A. Source Galaxies 


l. Shapes 


The DES Y3 shear catalog is derived using the META- 
CALIBRATION pipeline [50][151], which infers the ellipticity 
and similar photometric properties of objects using informa- 
tion from the ri z-bands. The pipeline is similar to that used 
in the DES Y1 analysis [152], but with a number of updates, 
including improved PSF solutions [118], improved astromet- 
ric solutions (115], and the inclusion of an inverse-variance 
weighting for the galaxies. METACALIBRATION is able to 
self-calibrate the initial estimate of the shear field from the 
measured galaxy shapes, including sample selection biases. 
The current METACALIBRATION implementation, however, 
does not correct for a shear-dependent detection bias 
that is coupled with object blending effects, which we find to 
cause a multiplicative bias in the shear at the level of 2-396. 
This residual bias is calibrated using image simulations [120]. 
Objects are included in the catalog if they pass a number of 
selection cuts designed to reduce potential systematic biases 
[117]. After additional footprint masking to match the lens 
catalogs, the final DES Y3 shear catalog yields 100 million 
galaxies covering an area of 4143 deg?, with a weighted effec- 
tive number density neg = 5.9 per arcmin? and corresponding 
shape noise oe = 0.26. 


The catalog has passed a variety of empirical tests [117], 
mostly aimed at identifying residual additive biases in the 
shear estimates. Systematic errors related to PSF modeling 
were shown to be negligible for the DES Y3 analysis, due to 
improved PSF modeling [118]. The B-mode signal was also 
shown to be consistent with zero. Other tests included the 
dependence of the shear estimates on galaxy and survey prop- 
erties. 


While shear calibration is typically viewed as separable 
from redshift inference, which is described in the following 
section, we also account for the first time for how blending 
correlates the ensemble shear calibration in each redshift bin 
with corrections to the effective shape of the n(z) of each 
of four redshift bins [120]. These corrections stem from a 
blending-detection bias, which biases both the ensemble av- 
erage shear (some blends will only be detected as separate 
objects depending on the shear) and redshift distribution in 
a potentially correlated way. One way to treat these effects 
coherently is to model the multiplicative shear calibration as 
a scaling of the total number density in each redshift bin, 
and fit these effects fully in redshift space. We would model 
then the shear calibration bias as e = fj^ n4(z)7\(z) 4- c. 
for additive shear bias c, observed ellipticity e, true ellip- 
ticity y, shear component j, and redshift bin 2. In prac- 
tice we continue to separate a scalar multiplicative (m) shear 
bias component to be compatible with existing codes, where 
n~(z) x (1+ m)n(z). 


MagLim n(z) cocoate 


RedMaGiC n(z) 


FIG. 1. The source (top), MagLim lens (middle), and redMaGiC 
lens (bottom) redshift distributions. The histograms are normalized 
to integrate to the total weighted galaxy density (arcmin ?) in each 
tomographic bin. The equivalent 1c uncertainties on the redshift dis- 
tributions are indicated by the shaded regions. The distributions have 
been corrected by non-zero mean and width offsets derived in the rel- 
evant photo-z uncertainty models. We adopt MagLim as our fiducial 
lens sample in this work, and use only redshift bins 1—4. 


2. Photometric redshifts 


The full redshift inference process for the source galaxies 
relies on connecting information about deep-field galax- 
les to those in the wide field that are used for the cosmological 
analysis [154][155]. All galaxies in the deep fields with similar 
properties are clustered together into different ‘phenotypes’ 
via a self-organizing map (SOM), while the same is also sep- 
arately done for all galaxies in the wide field. The deep-field 
galaxies have much lower photometric noise and additional 
wavelength information (i.e., overlapping infrared photome- 
try), so more phenotypes can be uniquely identified than for 
the wide-field galaxies. A redshift distribution is inferred for 
each of the deep-field galaxy phenotypes using overlapping 
spectroscopic [1561/1160], and photometric COSMOS 
and PAUS redshift measurements. We then create a 
probabilistic mapping between the deep- and wide-field phe- 
notypes using the Balrog simulation [121]. For example, if a 
given wide-field galaxy phenotype was mapped uniquely onto 


a single deep-field galaxy phenotype, its redshift distribution 
would be determined by the available redshift measurements 
of the deep-field galaxies that share that particular phenotype. 
In practice, the mapping is much more complicated: each 
wide-field phenotype has a non-zero probability of coming 
from many deep-field phenotypes, but the algorithm for gen- 
erating an n(z) for that galaxy phenotype is simply a weighted 
average. A given redshift bin is defined by a unique subset of 
many wide-field galaxy phenotypes, and its n(z) then follows 
by averaging over these phenotypes. The four source redshift 
bins have edges z € [0.0, 0.36, 0.63, 0.87, 2.0]. 

The process we use to account for uncertainties accumu- 
lated in each step of this process is summarized in Ref. [123]. 
These are due in part to shot-noise and cosmic variance in the 
redshift samples and deep fields [163], and photometric cali- 
bration uncertainty. At low redshift, this uncertainty is primar- 
ily due to uncertainties in the photometric calibration, while at 
high redshift it is due to a combination of cosmic variance and 
uncertainties in the redshift samples. Uncertainty in the n(z) 
due to these effects are modeled or measured, and we generate 
many realizations of the redshift distribution, n;(z), that ap- 
propriately sample the joint space of this uncertainty without 
relying on a simple parameterization like mean and width. 

The emerging set of redshift realizations suffer from one 
further source of uncertainty that has not been explicitly in- 
cluded before: blending. Galaxies that are nearby one an- 
other when projected on the sky can actually be very far apart. 
Detection and measurement algorithms can misinterpret these 
blends and report not only incorrect shapes but also incorrect 
number densities or redshifts. To account for this effect, we 
created realistic simulations and apply the same detec- 
tion and measurement pipeline used for the DES data to ob- 
tain the “observed” number density, shape, and photometric 
redshift of a known simulated object population that matches 
our deep field data, from which the impact of blending can be 
inferred. The result is a likelihood model describing the im- 
pact of blending on the joint shear calibration and n(z) shape 
that we add to the n;(z) [120]. 

We empirically constrain the likelihood of each n;(z) us- 
ing information from galaxy clustering on small scales that is 
not used in the primary 3x2pt observables [125]. We know 
that galaxies are likely to be found near other galaxies due 
to gravitational clustering. Therefore, if there is a galaxy in 
a given direction whose redshift is known, it is likely that 
nearby galaxies on the sky are at a similar redshift. We make 
use of multiple galaxy samples with well-determined redshifts 
and cross-correlate them with the wide-field source sample, 
thereby obtaining a likelihood for each of the n;(z), which is 
jointly sampled with the models that produce the n;(z) before 
we account for blending effects. We produce several thou- 
sand n;(z) samples. We show the final redshift distributions 
and their uncertainties in the top panel of Fig. 

To sample over n;(z) in a likelihood analysis, we introduce 
a set of hyper parameters to our model that rank the n;(z) 
in multiple dimensions [127]. Ref. has demonstrated 
that our constraints on the variation in n;(z) (i.e., shown in 
Fig.[I) are sufficiently precise that uncertainty in higher-order 
modes in the n;(z), besides the mean redshift, were not ex- 


pected to impact our cosmological constraints at a significant 
level in Y3. Thus in practice we simply sample the mean of 
the redshift distribution in the four redshift bins within a Gaus- 
sian prior based on the measured variance in the mean of each 

Finally, the measured two-point functions themselves fur- 
ther constrain the possible values of the redshift distributions 
via self-calibration in 3 x2pt. We further augment this by ex- 
plicitly using a set of the ratios of the galaxy—galaxy lensing 
signal on small scales between source redshift bins sharing 
the same lens bin [128], which contains information not used 
in the standard 3 x2pt analysis. These scales are too difficult 
to model robustly in full, but the ratios are to first order inde- 
pendent of cosmological model and depend primarily on the 
redshift distribution and intrinsic alignment parameters, and to 
a lesser degree on any redshift dependent bias in shear calibra- 
tion. This small-scale shear ratio likelihood is jointly sampled 
in the cosmological analyses. 


B. Lens Galaxies 


We have selected two galaxy populations (MagLim and 
redMaGiC) that serve as ‘lenses’ in galaxy-galaxy lensing 
measurements and for galaxy clustering measurements. The 
fiducial results presented in this work use the MagLim sample. 
We now describe the two lens samples. 


1. MagLim sample 


We have selected a magnitude-limited lens sample 
141], which results in 10.7 million galaxies. This “MagLim” 
sample is defined with a magnitude cut in the 7-band that de- 
pends linearly on redshift, ; < 4z + 18, where z is the pho- 
tometric redshift estimate from the Directional Neighborhood 
Fitting (DNF) algorithm [115][164]. This selection was opti- 
mized for wCDM constraints (135]. 


The MagLim sample is divided into six tomographic 
bins from z = 0.2 to z = 1.05, with bin edges z = 
[0.20, 0.40, 0.55, 0.70, 0.85, 0.95, 1.05]. The redshift distribu- 
tions from DNF are shown in the middle panel of Fig. [I]and 
have been validated using galaxy clustering cross-correlations 
[122]. Weights are derived to account for correlations in the 
number density with survey properties [136]. Further valida- 
tion and characterization of the sample is described in Refs. 
[141]. After unblinding, we discovered issues with the 
sample above z — 0.85, which lead to disagreement between 
the galaxy clustering and galaxy-galaxy lensing signal, and 
contribute to a substantially poor model fit to any cosmolog- 
ical models considered in this work (i.e., the two right-most 
panels of Fig. [2). This led us to remove these redshift bins in 
the fiducial analysis, which is discussed further in Secs. 
& 


2. redMaGiC sample 


This sample is selected with the redMaGiC algorithm 
[165], which results in 2.6 million galaxies. redMaGiC selects 
Luminous Red Galaxies (LRGs) according to the magnitude- 
color-redshift relation of red sequence galaxies, calibrated 
using spectroscopic redshifts. The sample has a luminosity 
threshold Lmin and approximately constant comoving den- 
sity. The redMaGiC sample has approximately 3096 narrower 
redshift distributions than MagLim, but approximately one- 
fourth the number of objects. 


We split the redMaGiC sample into five tomo- 
graphic bins, selected on the redMaGiC redshift 
point estimate quantity. The bin edges used are 
z = [0.15,0.35,0.50, 0.65,0.80,0.90]. The first three 


bins use a luminosity threshold of Lmin > 0.5L, (the ‘high 
density’ sample). The last two redshift bins use a luminosity 
threshold of Lmin > 1.0L. (the ‘high luminosity’ sample). 
The redshift distributions are computed by stacking samples 
from a non-Gaussian redshift PDF of each individual red- 
MaGiC galaxy. Each distribution is built from several draws 
of the redshift PDF and are shown in the bottom panel of 
Fig.[1] The mean and RMS width of the redshift distributions 
are validated using galaxy clustering cross-correlations in 
Ref. [122]. 

Weights are derived to account for correlations in the num- 
ber density with survey properties [136]. Further valida- 
tion and characterization of the sample is also described in 
Refs. [140]. We find a potential residual systematic in 
the redMaGiC sample at all redshifts, which does not impact 
ACDM inference and is also discussed in Sec.[V C] 


Ill. TWO-POINT MEASUREMENTS 


To extract cosmological information from the lens and 
source catalogs, we compute three sets of two-point corre- 
lation functions, which each measure information about how 
mass in the Universe is clustered. There are two fields rep- 
resenting the matter distribution that we can access with a 
galaxy survey: 1) the galaxy density field and 2) the weak 
lensing shear field. These two fields lead to these three sets of 
measured two-point functions. 

Galaxy Clustering: The two-point function between lens 
galaxy positions in redshift bins i and j, w'7(0), describes 
the excess (over random) number of galaxy pairs separated 
by an angular distance 0. The estimator for w^" (0) and its 
measurement and validation process are described in detail in 
Ref. [136]. We only use the auto-correlations of the measured 
w? (8) in our analysis; these are shown with their uncertainties 
in Fig. ]for MagLim and in App. [C]for redMaGiC. 

Galaxy—Galaxy Lensing: The two-point function between 
lens galaxy positions and source galaxy tangential shear in 
redshift bins i and j, 2d (0), describes the over-density of 
mass around galaxy positions. The matter correlated with the 
lens galaxy alters the path of the light emitted by the source 
galaxy, thereby distorting its shape. The estimator for y% (0) 
and its measurement and validation process are described in 
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FIG. 2. The measured w(6) correlation functions for each tomographic bin i of the MagLim lens galaxies (indicated by the 7,7 label in each 
panel). The best-fit ACDM model from the fiducial 3 x2pt analysis is plotted as the solid line in the top part of each panel, while the bottom 


part of each panel shows the fractional difference between the measurements and the model prediction, (ws 


—w")/ow (with y-axis range 
+50). In both the top and bottom part of each panel, 1c error bars are shown. Small angular scales where the linear galaxy bias assumption 


breaks down are not used in the cosmological analysis; these scales are indicated by grey shading. Bins 5 & 6 are not used in the final analysis. 
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FIG. 3. The measured +;(0) correlation functions for each tomographic bin combination using the MagLim sample. In each panel, the label 
i, j refers to MagLim lens tomographic bin 4 and the source bin j The best-fit ACDM model from the fiducial 3 x 2pt analysis is plotted as the 
solid line in the top part of each panel, with dotted curves indicating a negative model fit. The bottom part of each panel shows the fractional 
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of each panel, lo error bars are included. Small angular scales where the linear galaxy bias assumption breaks down are not used in the 


cosmological analysis; these scales are indicated by grey shading. Bins 5 & 6 are not used in the final analysis. 
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FIG. 4. The measured small-scale shear ratio values for each tomo- 
graphic bin combination using the MagLim sample, with lo error 
bars indicated. The x-axis identifies the two source bins that make 
up the measured ratio. The best-fit cosmological model from the 
fiducial 3 x 2pt analysis is over-plotted as the solid line for each set 
of lens-bin shear ratios. 


detail in Ref. [137]. The measured ^; (0) and their uncertain- 
ties are shown in Fig. 3]for MagLim and in App. [C]for red- 
MaGiC. In addition, we include small-scale shear ratio infor- 
mation below the scale cuts used for y+. These ratios are con- 
structed from ^, measurements using different source galaxy 
bins, while keeping the lens bin fixed. This effectively erases 
their dependence on the galaxy power spectrum, but keeps 
information about redshift calibration, shear calibration, and 
galaxy intrinsic alignment. A detailed description of shear ra- 
tios and their validation can be found in Ref. [128]. Fig. 
shows the shear-ratio measurement and uncertainties. This 
shear ratio data is included when analyzing all combinations 
of the three primary two-point functions in our analyses, un- 
less otherwise noted. 

Cosmic Shear: The correlation between source galaxy 
shears in redshift bins ? and j is described by the two functions 
£ i: (0), which are the sum and difference of the products of the 
tangential- and cross-components of the projected shear. The 
estimator for £e (0) and its S T and validation pro- 
cess are described in detail in Refs. (142][143]. The measured 
gu (0) and their uncertainties are fe in Fig. 

The total data vector includes measurements from five or 
six lens redshift bins and four source redshift bins, shown 
in Fig.|1} split into 20 logarithmic angular bins between 2.5 
and 250 arcmin, for a total of 1300 elements (not includ- 
ing shear-ratio). After bin pair removal for w(0), impos- 
ing a post-unblinding maximum lens redshift cut, and other 
scale cuts, 462 elements remain in the final 3 x2pt data vec- 
tor. The scale cut choices and their validation are described 
in Refs. [143], but are generally set to con- 
trol the impact of unmodeled non-linear effects (e.g., bary- 
onic effects on the matter power spectrum or higher-order 
galaxy bias) to better than 0.30 in the Qm—Ss plane in A- 
and wCDM. A choice of scale cuts that require biases meet 
our fiducial requirements in ACDM-only (ACDM-optimized) 
leaves 508 data points. All measurements are made using 
TreeCorr [166]. We find a total signal-to-noise S/N = 87 
for the 3x2pt data vector after fiducial scale cuts, where 
S/N = ania Emeders V EmodelC~ Emodel with covari- 
ance matrix C and best-fit model £544. This is a factor of 
2.1 improvement over the DES Year 1 3x2pt S/N. 


All of these measurements are related to the underlying 
clustering of matter in the Universe, but in different ways. 
The relationship between the galaxy density and the underly- 
ing matter density is complex and needs to be modeled 
with care. Alternately, the shape distortions depend more di- 
rectly on the intervening matter, but the measurements them- 
selves — especially of shapes and the redshift distribution — 
require greater care. Our work on this calibration is summa- 
rized in Appendix|[A] The advantage of using all of these mea- 
surements is that the systematic difficulties differ from one to 
another, but they all measure the same underlying matter field. 
Hence, by comparing the results from each set, we obtain a 
measure of consistency and additional ability to self-calibrate 
systematics, thereby giving confidence that we are correctly 
inferring information about the clustering of matter and the 
cosmological model. 


IV. ANALYSIS 


To infer parameters p from the measured two-point func- 
tions, we compare data organized in a “data vector" D, 


= {0 (0), 57 (8), £7 ()}, (1) 


to a theoretical model prediction organized in a vector T; 
of two-point correlation functions that are computed using the 
parameters p of a given model M, 


Tu(p) = (w(6.p),a? (0, p) £2(.p), — 


assuming a Gaussian likelihood, 


n Taye 
£(Dlp, M) x e 3L T») c (b-T«()) (3) 
Here C is the data covariance, which is obtained through an- 
alytic modeling as described and validated in Ref. [130]. 

We construct a posterior probability distribution for the pa- 
rameters p of the theoretical model given the data D as 


P(p|D, M) x £(D|p. M)P(p|M), (4) 


where P(p|M) is a prior probability distribution on the pa- 
rameters. The proportionality constant is given by the inverse 
of the Bayesian Evidence 


P(D|M) = I dp L(D|p,M)P(p|M), 6) 


which corresponds to the marginalized probability of a dataset 
being produced under a given theoretical model. 

This section summarizes the theoretical model and param- 
eterization we use for T m (p), which is described in more de- 
tail in Ref. and validated in Refs. [119] [129] [139] [140]. 
For clarity, we drop the parameter argument of the theoreti- 
cal model predictions, such that, e.g., the predicted clustering 
signal is simply denoted as w? (0). 

We report the mean in each parameter, along with the 68% 
confidence limit (CL) of posterior volume around the mean. 
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FIG. 5. The measured £4 (0) correlation functions for each tomographic bin combination, with labels as described in Fig.|3| The best-fit 
ACDM model from the fiducial 3 x2pt analysis is plotted as the solid line in the top part of each panel, while the bottom part of each panel 
shows the fractional difference between the measurements and the model prediction, (esp — £2)/oe, (with y-axis range +50). In both 
the top and bottom part of each panel, 1c error bars are included. The shaded regions (both light and dark) indicate scales not used in the 
fiducial analysis, primarily due to uncertainties in the impact of baryonic effects. The lighter shaded regions indicate scales that are used in an 
ACDM-optimized analysis, which meets our criterion for scale cuts described in Sec.[IV]in ACDM only. 


For completeness, we also report the best-fit maximum pos- 
terior values. We have used both a parameter-level and x? 
criterion for limiting the contribution of any systematic error 
to bias in the cosmological parameters. The threshold for this 
criterion is intended to limit the expected total bias in the 2D 
marginalized (2,,—5s plane from several independent poten- 
tial sources of model bias to be contained within the 68% C.L. 
region (< 0.30 for any single contribution). The differ- 
ence between the mean and best-fit values can give an indi- 
cation of the magnitude of projection or non-Gaussian effects 
in the marginalized parameter posteriors. The estimated im- 
pact of projection or volume effects in the DES Year 3 3x2pt 
posteriors are tested and summarized in Ref. (129]. We also 
provide a 2D figure of merit (FoM) defined for two parame- 
ters as FoM,, », = (det Cov(pi,po)) /? [167] [168]. The 
FoM is proportional to the inverse area of the confidence re- 
gion in the space of the two parameters, and can be considered 
a summary statistic that enables a straightforward comparison 
of constraining power of experiments or analysis scenarios. 


The analysis was designed and validated without access to 
the true cosmological results to protect against confirmation or 


Observer bias. This process is described in detail in App.[D] 


A. Model 


We model the observed projected (lens) galaxy density con- 
trast ô’ (ñ) as a combination of projected galaxy density 


contrast and modulation by magnification, ó PT 


ops (lt) = 5, (8) + ô, (ñ) (6) 
for position vector ñ, where i and j represent the redshift bin. 
The observed shear signal y is modeled as the sum of gravita- 
tional shear, ^c, and intrinsic alignments, ej, 

WÂ) = 6, (8) + Ea À), (7) 
with o the shear components. While B-modes produced by 
higher-order weak lensing effects are negligible for our analy- 
sis, it is important to account for B-modes generated by intrin- 
sic alignments in the computation of cosmic shear two-point 


TABLE I. The model parameters and their priors used in the fiducial 
flat ACDM and wCDM analyses. The parameter w is fixed to —1 in 
ACDM. The parameters are defined in Sec. 


Parameter Prior 
Cosmology 
Qu Flat (0.1, 0.9) 
10? A, Flat (0.5, 5.0) 
Ns Flat (0.87, 1.07) 
Qb Flat (0.03, 0.07) 
h Flat (0.55, 0.91) 
1090, 2? Flat (0.60, 6.44) 
w Flat (—2.0, —0.33) 
Lens Galaxy Bias 
bi(i € [1, 4]) Flat (0.8, 3.0) 
Lens magnification 
Cei Fixed 0.42 
c Fixed 0.30 
CP Fixed 1.76 
Cr Fixed 1.94 
Lens photo-z 
Az x 10? Gaussian (—0.9, 0.7) 
Az? x 10? Gaussian (—3.5, 1.1) 
Az} x 10? Gaussian (—0.5, 0.6) 
Az} x 10? Gaussian (—0.7, 0.6) 
cli Gaussian (0.98, 0.06) 
02 Gaussian (1.31, 0.09) 
o? Gaussian (0.87, 0.05) 
o Gaussian (0.92, 0.05) 
Intrinsic Alignment 
ai (i € [1,2]) Flat (—5,5) 
m Ge [1, 2]) Flat (—5,5) 
bra Flat (0, 2) 
Zo Fixed 0.62 
Source photo-z 
Azi x 10? Gaussian (0.0, 1.8) 
Az x 10? Gaussian (0.0, 1.5) 
Az? x 10? Gaussian (0.0, 1.1) 
Az x 10? Gaussian (0.0, 1.7) 
Shear calibration 
m! x 10? Gaussian (—0.6, 0.9) 
m? x 10? Gaussian (—2.0, 0.8) 
ue x 10? Gaussian (—2.4, 0.8) 
må x 10? Gaussian (—3.7, 0.8) 


correlation functions. 
can be written as 


In Fourier space, this decomposition 


Eds, AOO, ©) 


with the convergence field 


KI (fA) = dx W(x) 6m (ÂX, x) ; (9) 


where ôm is the 3D matter density contrast. The galaxy den- 
sity contrast ô is related to ôm via a linear galaxy bias b;. The 
tomographic lens efficiency is 


dU HS 4 !— 
ee i darax). (10) 
2 x a(x) 


wio) = 


11 


x is the comoving distance, x; the comoving distance to the 
horizon, n; (x) the source galaxy number density distribution, 
and a(x) the scale factor. 


I. Two-point statistics 


The angular power spectra C(£) of these observed fields 
can be written as 


a (6) —CE, (0) + C3 (0) + C4, (0 + Czy, (0) 
EM "es 
e (£) = o )+ Cj KORRO LU) P J ale ) 
(4) 


£) 2C 5, (6 )+ 2 Sn Qs oa (4) 


E CP. (£) + 20% esp Ad F 20$ snos, (l P 

(11) 

With the exception of the galaxy clustering power spectra 
C'5,,.5,,," Which are evaluated using the method described in 


Ref. [169], we calculate the angular cross-power spectrum be- 
tween two fields A, B using the Limber approximation 


£i 
2x) , 
X 
(12) 
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with PAp the corresponding three-dimensional power spec- 
trum, which is specified by the parameterization choices sum- 
marized in The kernels WA B correspond to W3 for 


shear and the lens galaxy density n; for position. The two- 
point correlation functions within an angular bin [Onin, Omax] 
are related to the projected power spectra as 


)= 2,90 (6 
= Tat 
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Omin, Omax) OF done (4) 


Omin, Omax) E qu (13) 


dobs 


with Gn analytic functions detailed in Refs. |129] . 


B. Parameterization and Priors 


We sample the posterior of these measurements in two cos- 
mological models: flat ACDM and wCDM, with the sum of 
the three neutrino masses as a free parameter, where the im- 
pact of neutrino mass on the power spectrum is modeled via 
a fitting function [? ]. ACDM contains three energy densities 
in units of the critical density: the total matter density Qm, 
the baryonic density Q», and the massive neutrino density 2. 
We vary Q,h?, where h is the Hubble parameter, as a free 
parameter, while noting that it is often fixed in other cosmo- 
logical analyses to be zero or to the minimum mass allowed 
by oscillation experiments m, = 0.06 eV [170]. 


The other cosmological parameters we vary within ACDM 
are the Hubble parameter h, the amplitude of primordial scalar 
density perturbations A,, and the spectral index n, of the 
power spectrum. We assume a flat model, with Q, = 1— Qm. 
In wCDM, we allow for a free dark energy equation-of-state 
parameter w that is constant in time (in ACDM, this is fixed 
to w = —1, corresponding to a cosmological constant). Thus 
ACDM includes six free cosmological parameters and wCDM 
contains seven. The prior ranges for cosmological parameters 
in Table|[IJare either motivated by physical constraints (e.g., an 
accelerating universe requires w « —1/3), or for parameters 
that are not strongly constrained by the DES data, typically 
given a range that encompasses five times the 68% C.L. from 
relevant external constraints. In analyses that sample exter- 
nal CMB likelihoods, we include the optical depth 7 as a free 
parameter. 

We will typically refer to the amplitude of density pertur- 
bations at z = 0 in terms of the RMS amplitude of mass on 
scales of 8h~! Mpc in linear theory, ca. The constraints on 
the amplitude and density of matter fluctuations are degener- 
ate in our analysis, and we will also refer to the parameter 
Sg, which describes the width of the posterior in the direction 
roughly orthogonal to the primary degeneracy direction for 
cosmic shear in the ca-(,, plane [171], though this does not 
hold exactly for 3x2pt and changes with effective redshift. 

In addition to these cosmological parameters, our fiducial 
analysis includes an additional 25 free parameters, for a to- 
tal of 31 (32) parameters in ACDM (wCDM). These addi- 
tional parameters describe astrophysical and systematic con- 
tributions to the measured signal. The effective linear galaxy 
bias of lens galaxies in each redshift bin is parameterized 
by a scalar b/. We also test and apply a nonlinear galaxy 
bias model (with one extra free parameter per redshift bin) 
[141], which is described in App. The intrin- 
sic alignment of galaxies is modeled with the Tidal 
Alignment and Tidal Torquing (TATT) model (174], which is 
parameterized by an amplitude a; and redshift power-law n; 
parameter (with redshift pivot zọ = 0.62), for each of the (1) 
tidal alignment- and (2) tidal torquing-sourced terms in the 
model, as well as an effective source galaxy bias parameter 
bra, which is described in further detail in Refs. [129] [174]. 
The TATT model contains the commonly employed non- 
linear linear alignment (NLA) model in the ag = bra = 0 
subspace. The amplitude of the lens magnification term in 
Eq. [6] depends on the slope of the lens sample’s luminosity 
and size distribution at the sample detection limit. The cor- 
responding parameter CÝ is calibrated from the data, as de- 
scribed in Ref. [139], and held fixed to that value. Nonlocal ef- 
fects in ^ can significantly contaminate larger angular scales 
with nonlinear information due to integration of the projected 
mass within a given angular separation from the center of the 
halo. This is mitigated by analytically marginalizing over a 
free point-mass contribution to ^; in all analyses [175]. 

Photometric redshift systematics are parameterized by an 
additive shift to the mean redshift of each bin, Az for lenses 
and Az? for sources, where the true redshift distribution is 
related to the photometric redshift distribution np, such that 


ni(z) = DA C — As). (14) 
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In addition, differences in the width of the lens redshift distri- 
bution are important at DES Y3 precision, which we parame- 
terize by a stretch o}, such that 


n! (z) = oin (ot [z — (z)] + (2)) : (15) 


Finally, uncertainty in the shear calibration bias is parameter- 
ized by m^, where the measured ellipticity e; is related to the 
true shear yj in each bin by 


ej = (1+ m')y. (16) 


The full set of parameters (cosmological, astrophysical, and 
systematic) and their priors are summarized in Table [i] 
Differences in the redMaGiC analysis are described in 


App. 


C. Likelihood Analysis 


Our likelihood analysis uses two independently developed 
analysis and inference pipelines, COSMOSIS and Cos- 
MOLIKE [177], which have been validated against one another 
to ensure they produce consistent predictions of the observ- 
ables and final cosmological constraints. A comparison of the 
theory predictions from COSMOSIS and COSMOLIKE is pre- 
sented in Ref. [129]. The residual offset of x? < 0.2 between 
3x2pt model data vectors obtained from both codes in this 
analysis at a reference cosmology is found to have negligible 
impact on parameter constraints and we conclude that both 
pipelines can be used interchangeably. 


1. CosmoSIS 


This pipeline uses the CAMB Boltzmann code 
to compute underlying background quantities and the linear 
matter power spectrum, and the HALOFIT version pre- 
sented in Ref. for the non-linear power spectrum. It 
then generates theory predictions following the model de- 
scribed in section [IV A] and using the Fast-PT method 
for non-linear galaxy bias and the TATT model for intrinsic 
alignments. Non-Limber integrals are computed following 
the method of [169]. Accuracy parameters throughout the 
pipeline are chosen by requiring the log-likelihood to differ 
by less than 0.05 from a high precision calculation. For chains 
including Planck CMB measurements [13], we use the Planck 
2018 public likelihood code [183] 

The version of CosMOSIS?| used for the analysis may be 
found in the des-y3 branch of the repositories. The Cos- 
MOSIS runs presented here use the PolyChord sampling 


5 While running the Planck CLIK lensing likelihood from plc-3.0, discrep- 
ancies between the constraints obtained using the likelihood code and the 
publicly released chains were found. The disagreement has been identified 


to originate from the treatment of the linear correction term to the theory 
oo 


Cy” spectrum. This has been corrected in the upstream plc-3.01 release. 


method [185], for both posterior samples and Bayesian 
evidence. The shear calibration values m; are used as fast 
parameters. The PolyChord parameters we use for our 
fiducial runs are: fast fraction = 0.1, live points = 500, 
num repeats = 60, tolerance = 0.01, and boost. posteriors 
— 10.0. 

Analyses with the COSMOSIS pipeline also use the non- 
Gaussian covariance matrix from COSMOLIKE described be- 
low. 


2. CosmoLike 


This pipeline uses the CLASS Boltzmann code [186] to 
compute underlying background quantities and linear and 
non-linear matter power spectra, using the HALOFIT version 
presented in Ref. for the latter. The theory predictions 
are calculated using the model described in section [IV A] re- 
lying on the FAST-PT method to evaluate integrals 
over perturbation theory kernels for non-linear galaxy bias 
and the TATT intrinsic alignment model. The computation 
of non-Limber integrals for galaxy clustering further employs 
the FFTLog implementation of Ref. [169]. 

The evaluation time of angular two-point statistics in 
COSMOLIKE is optimized through a series of interpolation 
schemes, for which runtime-optimized accuracy settings are 
validated through comparison to high-accuracy evaluations 
with slow runtime. Most DES Y3 likelihood analyses with 
COSMOLIKE employ the EMCEE sampler, c.f. for 
a detailed comparison of sampler configurations for COSMO- 
LIKE likelihood analyses. 

The CosMoCov module of COSMOLIKE is used 
to generate covariances for DES Y3 analyses, which include 
Gaussian and non-Gaussian terms and account for the 
effect of the survey geometry on shape and shot-noise terms 
[101]. 


D. Tests on Simulations 


Our model and many other components of our analysis 
have been validated end-to-end on a suite of 18 cosmological 
simulationg] [1914193]. Validation is performed on the 
mean of measurements from all 18 of these simulations with- 
out shape noise, including photo-zs and marginalizing over 
all cosmological and nuisance parameters. We have verified 
that we can recover the correct cosmology with our fiducial 
analysis to within approximately 0.3c in the 2D og—OQm plane 
(ACDM) and w—,, plane (wCDM). We have shown that for 
a more stringent test in the absence of photometric redshift 
and shear calibration uncertainties (using true redshifts) our 
model is able to reproduce the mean £+ and w + y+ measure- 
ments from our 18 simulations with a x? of 1.4 for cosmic 


7 Each simulation assumes a flat ACDM cosmology with Qm = 0.286, 
Q» = 0.046, h = 0.7, ns = 0.96, and og = 0.82. 
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shear (207 data points), 4.5 for w(@) (53 data points), and 9.1 
for y+(0) (232 data points). These x? numbers are relative to 
the fiducial covariance for a single DES Y3 realizations, but 
with a measurement that is the average of 18 realizations with- 
out shape noise. Thus, they represent the potential systematic 
X? contribution due to model inaccuracies, and shouldn't be 
interpreted as a goodness-of-fit metric. 

In addition to these model tests, we have also investigated 
the systematic uncertainty inherent to our redshift inference 
process using these simulations. We have shown that the 
three independent source redshift n(z) estimates - SOMPZ, 
source-lens clustering, and shear ratios — produce consistent 
constraints on the redshift distribution. We have also per- 
formed our fiducial analysis using a redMaGiC-like lens sam- 
ple, assuming source redshift distributions that are calibrated 
using the same three methods and lens redshift distributions 
estimated from redMaGiC. We found that the constraints from 
this analysis are consistent with those that use the true red- 
shift distributions from the simulation. The final constraining 
power is similar between redMaGiC and MagLim, so we do 
not repeat the simulated analysis twice. 


E. Quantifying internal and external consistency 


To quantify consistency of internal and external data sets, 
we define a priori a process to guide decisions and conclu- 
sions before seeing the cosmological constraints. For internal 
consistency, we calculate the Posterior Predictive Distribution 
(PPD) and derive a (calibrated) probability-to-exceed p. 
In short, the idea is to draw realizations of a particular subset 
of the data vector for model parameters drawn from the pos- 
terior of the same subset (goodness-of-fit tests) or a disjoint 
subset (consistency tests). These realizations are then com- 
pared to actual observations and a distance metric is computed 
in data space, which is then used to compute the p-value. We 
test the goodness-of-fit for the two combinations of two-point 
functions, £x and w + y+, and the combination of all three 
after confirming they are mutually consistent. In all cases, we 
require as part of the unblinding criteria defined in App.[D]that 
p > 0.01. The validation of the use of PPD for these tests is 
described in Ref. [132]. 

To quantify consistency with external experiments, we have 
explored a variety of metrics in order to calibrate expectations, 
which are described in Ref. [133]. We studied the particu- 
lar case of quantifying consistency between DES 3x2pt and 
Planck CMB using both simulated DES Y3 data and real Y1 
data. The consistency metrics can be divided into two cate- 
gories: parameter-based, which measure relative deviations in 
the multi-dimensional parameter space, and Evidence-based, 
which also account for how well the individual and combined 
data sets fit the model. We discuss results in terms of at least 
one metric from each category: the parameter difference and 
Suspiciousness [194], along with the Evidence ratio to com- 
pare to DES Y1. These metrics can produce a probability-to- 
exceed, and we require the same criterion of p « 0.01 to con- 
clude there exists evidence for inconsistency between probes. 

Detailed results from these consistency studies are shown 
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FIG. 6. Marginalized constraints on the three parameters og, Sg = 
C8 V/ Nm /0.3, and Qm in the ACDM model from cosmic shear (£4 
blue), galaxy clustering and galaxy-galaxy lensing (y+ + w(0), or- 
ange) and their combination (3x 2pt, solid black). We also show 
a ACDM-optimized 3x 2pt analysis that is valid for ACDM using 
smaller angular scales in cosmic shear (dashed black). The marginal- 
ized contours in this and further figures below show the 68% and 
95% confidence levels. The top and side panels show 1D marginal- 
ized constraints with the 68% confidence region indicated. 


in App.[F] 


V. DES Y3 RESULTS: PARAMETER CONSTRAINTS 
A. ACDM 


The principal cosmological test of the DES Y3 data is to 
compare our data to the currently favored ACDM model. The 
model has six cosmological parameters, but also 25 nuisance 
parameters for a total of 31 free parameters (listed in Table |I) 
in our fiducial analysis with the MagLim lens sample. Recall 
also that, in the fiducial analysis, we use all four source-galaxy 
bins, but only the first four (out of six total) MagLim lens- 
galaxy bins. 

We first concentrate on comparing two subsets of 3x2pt 
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measurements: those from cosmic shear (£+) and those from 
the combination of galaxy clustering and galaxy-galaxy lens- 
ing (w + ^g). It is logical to compare these two subsets be- 
cause each can constrain the ACDM parameters to a similar 
precision, yet the constraints from the two subsets contain in- 
dependent information, since they are sensitive to the under- 
lying matter density fluctuations in different ways. Moreover, 
comparing the two subsets of the 3x2pt measurements pro- 
vides an internal consistency test. 

We must first check that the £+ and w + ^»; measurements 
are a good fit to the data, mutually consistent, and that their 
combination (3x2) is a good fit to the model. For each of 
the two measurements individually, we find a PPD result for 
model goodness-of-fit p(£+) = 0.21 and p(w + y+) = 0.02. 
The PPD result for consistency between the two model con- 
straints is p(£4 |w 4- y+) = 0.30, meaning that it is appropriate 
to combine £4. and w + y+. The joint 3x2pt goodness-of-fit 
is p(&4 + y + w) = 0.04. Finally, the shear-ratio data has 
goodness-of-fit p = 0.03 in this joint best-fit model. All of 
these p values meet our original criterion of p > 0.01 defined 
in App.[D] 

The marginalized constraints from each probe and the 
3 x2pt combination in the parameter space spanned by ag, Sg, 
and Qm are shown in Fig.|6| This is also summarized in Table 
and Fig. [7| which show the numerical constraints on these 
three parameters. The DES Y3 3x2pt constraints on the key 
parameters are 


Ss = 776 hie (0.776) 


Qm = 0.339+6:631 (0.372) (17) 


Gg = 0.733+0:039 (0.696). 


The 3x2pt contours in these parameters are not centered on 
the overlap of £+ and ^; + w due to degeneracies in the higher 
dimensional parameter space. 

We also perform two alternative 3x2pt analyses that use 
smaller scales. First, we perform a ACDM-optimized analysis 
that includes smaller-scale information in cosmic shear. This 
analysis meets our parameter bias requirements in ACDM 
(i.e., Sec.[IV]. but not wCDM. The 3 x 2pt results from the op- 
timized analysis are shown in the row labeled ‘ACDM-Opt.’ 
in Table |II| and Fig. The optimized results are consistent 
with the fiducial analysis, but are about 3096 more constrain- 
ing in the 2D marginalized 22,,—cg plane. The second alterna- 
tive analysis utilizes a more complicated nonlinear bias model 
in order to model smaller scale information in ^; + w(0), and 
is described in App. |E| The 3x2pt results from the nonlinear 
analysis, shown in the row labeled ‘NL bias’ in Table |I| and 
Fig.|7| are consistent with the fiducial analysis and lead to an 
increase of 15% in constraining power in the Qm—øs plane. 

While we found no significant evidence of internal incon- 
sistency with the ACDM model using the final MagLim lens 
selection, we have identified potential systematic modes in 
the data at high redshift for the MagLim sample and at all 
redshifts for the redMaGiC sample. We had agreed before 
seeing any cosmological results that we would pursue poten- 
tial systematics in the case where the results failed to suffi- 
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FIG. 7. Summary of marginalized constraints (mean and 68% CL) and maximum posterior values (crosses) on Sg, Qm, and og in ACDM. 
‘Ext. Low-z’ data consists of external SNe Ia, BAO, and RSD, while ‘All Ext.’ data consists of external SNe Ia, BAO, RSD, and Planck 
CMB with lensing. The top section shows constraints using only DES data, the middle section only external data, and the bottom section 


combinations of DES and external data. 
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FIG. 8. Constraints on the galaxy bias (b,) and effective intrin- 
sic alignment (IA) amplitude from tidal alignment (a1) and tidal 
torquing (a2) are shown per redshift bin. Constraints using both lens 
samples (MagLim and redMaGiC) are shown. The galaxy bias is 
expected to be different for both lens samples, but the IA amplitude 
constraints, which are a property of the source galaxy sample, are 
consistent. We do not necessarily expect a1 and a» to be consistent 
with one another. We sample over a power-law evolution of the IA 
amplitude, so the redshift evolution is forced to be smooth in a;. 


ciently fit any of the models considered in this work (ACDM 
and wCDM) at p « 0.01. Including MagLim lens bins 5 
and 6 caused a very poor model fit to both models, with 
p £: 5x 10-7. Based on this criterion, we applied a high-z cut 
to limit the MagLim sample to approximately the same red- 
shift range of redMaGiC post-unblinding. This change is dis- 
cussed further in App. [D] The two lens samples are compared 
and further details of this are discussed in Sec. but all 
issues that have been uncovered appear to be mostly orthog- 
onal to the 3x2pt ACDM parameter dimensions — that is, 
they do not significantly impact the inferred cosmological pa- 
rameters, and the cosmological parameters inferred from the 
two lens samples are consistent. This resilience of the 3x2pt 
combination of data and its ability to self-calibrate potential 
systematics in a subset of the two-point functions is one of 
the main motivations for pursuing this cosmological probe for 
large-scale structure. 

We find that the DES Y3 3x2pt analysis is able to add in- 
formation beyond the prior for 15 parameter dimensions in 
the model, three of which are cosmological. The cosmologi- 
cal modes that DES 3x2pt most improves with respect to the 
prior are obtained with the Karhunen-Loéve decomposition of 
the posterior and prior covariance, and are: 


— 0.77 _ +0.015 
pi = osha = 0.317 ^ 6014; 


p2 = maz 049 71 (18) 


ps = hnl 4ng 39 = 2,114945, 


The combined 3 x 2pt data is also able to simultaneously con- 
strain a variety of ‘astrophysical’ parameters that encode how 
galaxies are connected to the underlying dark matter perturba- 
tion field, namely the linear and nonlinear bias parameters and 
intrinsic alignment of galaxies. Constraints for these model 
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FIG. 9. Marginalized constraints on the two parameters Qm and w 
in the wCDM model from DES Y3 3x2pt. A dotted line indicates 
w = —] as given by the cosmological constant. 


parameters are shown in Fig. [8] We find slightly higher galaxy 
bias constraints for redMaGiC galaxies than in the DES Y1 
analysis using a similar redMaGiC sample. We find a prefer- 
ence for a slightly smaller intrinsic alignment amplitude than 
DES Y1. This value is consistent with the DES Y1 analysis, 
but is also consistent with zero intrinsic alignment. 


B. wCDM 


We also fit our data to the wCDM model, in order to test for 
evidence that the dark energy equation of state departs from 
its cosmological-constant value of w — —1. In wCDM, the 
dark energy density evolves with time with a constant w, such 
that ppp « (1 + z)9 (0), We show marginalized parameter 
posteriors for this model in Fig. p] and parameter values in 
Table []and Fig. 

We find similar levels of agreement between £4 or y;+w(0) 
as in ACDM, and a similarly good fit to the data within the 
wCDM model, but do not show constraints from these sub- 
sets of the data due to increased prior influence and parameter 
volume effects. The DES Y3 3x2pt constraint on the matter 
density and dark energy equation of state parameter are 


Om = 0.352+9:937 (0.339), 
w = —0.98+9:32 (—1.03). 


(19) 


To determine if there is a preference for the wCDM model 
over the ACDM model, we compute the Bayes factor 
P(D|ACDM 
p P(DIACDM) 


Ei und eine M Q0) 
P(D|wCDM) 
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A value of R greater than unity implies that the wCDM model 
is not favored. We find R = 4.3. This indicates that the late- 
universe large-scale structure probed by DES does not show 
evidence of needing the more complex dark energy density 
scenario of the wCDM model. We discuss further in Sec. [VII] 
more stringent tests of the ACDM model that leverage data 
across the age of the Universe. 


C. Lens sample comparison 


As optical surveys cover larger fractions of the sky and 
probe higher redshifts, photometric galaxy clustering be- 
comes both more powerful and more difficult to calibrate. Pre- 
vious DES analyses used a luminous red sample of galaxies 
with constant comoving density, redMaGiC. To ensure robust- 
ness, we pursued two lens samples for DES Y3: a magnitude- 
limited lens sample, MagLim, and the redMaGiC sample. The 
redMaGiC sample was optimized for better understood and 
smaller photometric redshift errors. The MagLim sample was 
optimized for wCDM constraints, balancing increased num- 
ber density vs. less well-constrained photo-zs, while allow- 
ing selection to higher redshifts than possible with redMaGiC. 
Comparing the inferred cosmological parameters of our mod- 
els from these two very different samples, which have fewer 
than 20% overlapping objects, allows us to infer potential un- 
corrected systematics from lens sample selection or photo-z 
calibration of the lenses. 

Measurements based on the second redshift sample, red- 
MaGiC, also have an acceptable overall model fit to the 
ACDM and wCDM models. The cosmic shear in this model is 
also consistent with the combination of galaxy-galaxy lensing 
and galaxy clustering. These measurements and model fits are 
shown in App. |C| We find a PPD result for model goodness- 
of-fit p(€+) = 0.25 and p(w + y) = 0.04, while the PPD 
result for consistency between the two model constraints is 
plé+|w + y) = 0.02. The joint 3x2pt goodness-of-fit is 
p(t + y + w) = 0.02. 

The marginalized constraints of each individual probe and 
the 3 x2pt combination on Sg and Qm in ACDM are shown in 
Fig.|11| while they are shown for Sg, w, and Qm in wCDM in 
an Both figures compare redMaGiC results to the fidu- 
cial 3x2pt using the MagLim sample. As described above, 
the cosmic shear or ^; + w(0) data alone are consistent with 
the 3x2pt model fit, though the ^; + w(@) data on their own 
prefer a smaller Sg value. This arises from the strong degen- 
eracy between og and galaxy bias in y + w(0). Alone, it 
prefers a lower value of og and higher value of galaxy bias. 
Adding cosmic shear information effectively fixes the value 
of Sg along that degeneracy, which brings the galaxy bias in 
3x2pt back down to a value more consistent with the DES 
Year 1 redMaGiC galaxy bias constraints. 

These results are consistent with the MagLim results and 
passed our unblinding requirements, including having a suffi- 
ciently good model fit to ACDM. However, after unblinding 
the results with redMaGiC we found evidence of internal ten- 
sion in the data. Because w(0) is not able to constrain cosmol- 
ogy on its own, this has limited impact on the combination 
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TABLE II. Summary of marginalized parameter constraints in ACDM. The mean and 68% CL are provided for each cosmological parameter, 
followed by the maximum posterior value in parentheses, except for neutrino mass, for which the 95% upper bound is given. Parameters that 
are not significantly constrained are indicated by a dash. All data have been re-analyzed with model and prior choices matching the DES Y3 
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FIG. 10. Summary of marginalized constraints (mean and 68% CL) and maximum posterior values (crosses) on Ss, Qm, and w in wCDM. 
‘Ext. Low-z’ data consists of external SNe Ia, BAO, and RSD, while ‘All Ext? data consists of external SNe Ia, BAO, RSD, and Planck 
CMB with lensing. The top section shows constraints using only DES data, the middle section only external data, and the bottom section 
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TABLE II. Summary of marginalized parameter constraints in wCDM. The mean and 68% CL are provided for each cosmological parameter, 
followed by the maximum posterior value in parentheses, except for neutrino mass, for which the 95% upper bound is given. Parameters that 
are not significantly constrained are indicated by a dash. All data have been re-analyzed with model and prior choices matching the DES Y3 
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FIG. 11. A comparison of the marginalized ACDM constraints of the 
two lens samples. Dashed contours show the cosmic shear (blue), 
galaxy-galaxy lensing and clustering (orange), and 3x2pt (black) 
constraints based on the redMaGiC lens sample. The 3x2pt red- 
MaGiC constraints marginalizing over a free Xiens parameter are 
also shown (dotted black), and the 3x 2pt MagLim constraints (solid 
black). The inferred cosmological parameters from 3 x 2pt are con- 
sistent in all three cases. 


of galaxy clustering and galaxy—galaxy lensing and no dis- 
cernible impact on the 3 x2pt combination in ACDM. Never- 
theless, it is important to understand the source of this internal 
tension in redMaGiC results and judge its impact on cosmo- 
logical inference. To do so, we modeled this inconsistency 
of the redMaGiC clustering and galaxy—galaxy lensing ampli- 
tudes with a systematic parameter Xjens, which is related to 
the connection of the galaxy—galaxy lensing and galaxy clus- 
tering two-point functions to the matter two-point function: 


y? (0) = POUR S n (0) 


where b; is the galaxy bias connecting the observable 7 or 
w(0) to the matter correlation function (Emm) or spectrum and 
X jens is the same for all redshift bins 7. We expect Xiens = 1 
in ACDM, if there are no systematic contributions to the sig- 
nals. The fiducial model described in earlier sections is thus 
identical to the model including in Eq. [21] with an additional 
constraint Xjens = 1. 

We show the result of marginalizing over a free Xjen;s in 
the redMaGiC 3x2pt analyses in Figs. [11]and [I2] We find 
a negligible impact on the primary cosmological parameters 
in ACDM, particularly Sg. We find Xjens = 0.87750 019. 
strongly inconsistent with Xi,,, = 1 in ACDM. If we fix 
Xiens to this value in the redMaGiC ^; + w(@) analysis, the 
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FIG. 12. A comparison of the marginalized wCDM constraints of 
the two lens samples. Dashed black contours show the 3x2pt con- 
straints based on the redMaGiC lens sample. The 3x2pt redMaGiC 
constraints marginalizing over a free Xjens parameter (dotted black) 
and the 3 x 2pt MagLim constraints (solid black) are also shown. The 
inferred cosmological parameters from 3x2pt are generally consis- 
tent, but in particular the redMaGiC results are sensitive to the impact 
of Xjens in w CDM, showing substantial shifts in the inferred param- 
eter values. 


contour in Fig. [LT shifts upward to agree with cosmic shear in 
Sg. The value of Xjens is correlated with the equation-of-state 
parameter w, so the redMaGiC wCDM constraint is strongly 
affected by this potential systematic. Adding the single free 
parameter Xj,,, in ACDM leads to an improvement in x? of 
25, while adding a free w leads to an improvement in x? of 7. 
Thus, Xjens clearly leads to a better model fit. 

After unblinding the results with redMaGiC, but before un- 
blinding those with MagLim, we decided to use the MagLim 
sample for our fiducial cosmological analysis if it showed no 
indication of this scale- and redshift-independent effect that 
is present in redMaGiC. This potential systematic was stud- 
ied at length between the unblinding of the redMaGiC sample 
and the MagLim sample. Studies of this effect are discussed 
in much more detail in (136][139][140]. We have demonstrated 
that the effect (and its relative impact vs. the clustering am- 


plitude of the MagLim sample) is roughly independent of red- 
shift, angular scale, or position in the survey footprint. 

After initial submission of this paper, we found that relax- 
ing the goodness-of-fit requirement for the red galaxy model 
selection in redMaGiC leads to a cosmological model fit con- 
sistent with Xj,,, = 1 and no significant change to the cos- 
mological parameter results. This test suggests that a color- 
dependent photometric issue is the source of Xjens! = 1, and 
is plausibly connected to background subtraction. A specific 
fix for this systematic at the image level has not been identi- 
fied, but these results pinning down the likely source of Xjens 
are described further in Ref. [140]. Further study of this ef- 
fect and pipeline modifications will continue for the final DES 
Year 6 analyses. 


1. Summary of possible non-systematic causes of Xiens # 1 


There are several classes of non-systematic explanations for 
X ens, all of which we believe are implausible given our data. 
These possible explanations are: 

Stochastic Bias: While the effect of Xjens on clustering 
and galaxy-galaxy lensing looks very similar to stochasticity, 
a decorrelation between the galaxy and matter distributions, 
predictions from galaxy bias models make this interpretation 
unlikely. In configuration space, perturbative stochastic terms 
are expected to contribute only at small separations r ~ R, 
(the Lagrangian size of halos), and to statistics that involve 
zero-lag correlators [67]. 

Lensing-is-low: Ref. reported that the 7, signal around 
luminous red galaxies is lower than expected from a model 
conditioned on their autocorrelation, which resembles the pat- 
tern seen in our redMaGiC sample. But with the possible 
exception of the large scale results in Ref. [196], the 
lensing-is-low result applies to models that fit to scales 
sensitive to complexities of the small scale dark matter-galaxy 
connection. There is still debate within the lensing-is-low lit- 
erature as to whether the effect can be accounted for by addi- 
tional complexity in these small scale models [197]. 

The DES Y1 results (which also used a redMaGiC sample) 
do not support the lensing-is-low scenario, nor do the DES 
Y3 results for MagLim in the redshift range of the redMaGiC 
sample. The DES Y3 results for the redMaGiC sample show 
what could be interpreted as galaxy-galaxy lensing being 10- 
15% percent lower than galaxy clustering at fixed cosmology 
(Planck 2015 in the case of Ref. [90]; DES 3 x 2pt cosmology 
in the DES Y3 results for Xjens). However, the more plausible 
cause is that the clustering of the Y3 redMaGiC sample is 
anomalously high, as indicated by internal consistency tests 
of the individual data vectors. While we are still studying 
the Xjens « 1 anomaly, we currently do not believe that it 
supports a conclusion that galaxy—galaxy lensing is “low.” 

Fundamental physics: Any dynamical modifications to ei- 
ther the Poisson equation or the shear equation generally 
changes the galaxy and matter distributions but their corre- 
lation is maintained, i.e. jens = 1 is maintained at linear 
scales. Beyond this possibility, any separation of the impact 
of relative “bias” between the two types of matter (apparent in 
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lensing vs. clustering) at the level of 15% would require sig- 
nificant fluctuations in the dark matter field, which would have 
substantial ramifications in other observables that we have not 
seen. Therefore, we conclude that a fundamental-physics ex- 
planation for Xjens < 1 would probably have to be very fine- 
tuned. 


2. Potential systematics in w(0) and Yt vs. Xiens Æ 1 


We now continue discussing the Xjens # 1 anomaly by 
comparing redMaGiC to MagLim, and commenting on poten- 
tial systematics in galaxy clustering and galaxy-galaxy lens- 
ing as the cause of the anomaly. 

We find the redMaGiC sample shows Xjens; < 1 at high 
significance at all scales and redshifts. The highest two red- 
shift bins of the MagLim sample, which have been removed 
from the analysis, also indicate Xjens < 1 at high significance, 
which is clearly visible in the model fit in those two bins of 
Fig.|2| In the redshift range overlapping the redMaGiC sam- 
ple, we find no evidence of a non-unity Xjens for MagLim. 
We discard the two high redshift bins for the MagLim sam- 
ple as a conservative choice. Based on our investigations so 
far and current understanding of theoretical extensions beyond 
wCDM, we do not believe these anomalies are indications of 
new physics. We have found plausible but unverified indica- 
tions that the origin may lie in potential systematics, e.g., as- 
sociated with the photometric uncertainty or background sub- 
traction for large or faint objects, or in the de-reddening pro- 
cess. 

These issues are the subject of ongoing investigations, 
which will be crucial for understanding photometric cluster- 
ing and its combination with galaxy-galaxy lensing in DES 
Y6 and beyond. However, while these measurements are po- 
tentially impacted at a level we can measure by some as yet 
unidentified systematic, this does not have a significant impact 
on ACDM cosmology when the three two-point functions are 
combined within 3 x2pt. For the MagLim sample our tests in- 
dicate that both the ACDM and wCDM constraints are robust. 
This self-calibration effect is one of the primary motivations 
for combining these different probes of the same underlying 
matter density field into the 3x2pt observable. 


VI. COMPARISON WITH OTHER DES DATA 


DES has produced competitive cosmological constraints 
using its four primary probes: galaxy clustering and weak 
gravitational lensing (3x2pt), type Ia supernovae (SNe Ia), 
galaxy cluster counts and masses, and BAO. Together, these 
probes have been demonstrated to provide dark energy con- 
straints that can be competitive with the best combined exter- 
nal constraints [107]. We describe each of them briefly below. 

Type Ia supernovae: The DES SNe Ia sample has 207 spec- 
troscopically confirmed SNe in the redshift range 0.07 < 
z < 0.85. The sample-building and analysis pipelines are de- 
scribed in a series of papers that detail the SN search and dis- 
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FIG. 13. A comparison of the marginalized constraints on parame- 
ters in the ACDM model from a variety of DES probes: large-scale 
structure and weak lensing (3 x 2pt; Y3 — black solid, Y1 reanalyzed 
— black dashed), type Ia supernovae (purple), galaxy cluster number 
counts and masses (orange and green), and BAO. The combination 
of DES Y3 3x2pt, SNe Ia, and BAO is shown in blue. Going from 
Y1 to Y3, we find approximately a factor of two improvement in the 
3x 2pt constraint in Qm-Ss plane. 


calibration [145] [202]; spectroscopic follow-up [203]; and se- 
lection bias [2041206]. The methodology and systematic un- 
certainties are found in Ref. [207]. These were used to con- 
strain cosmology and the Hubble constant [209]. These 
analyses included additional external low-redshift SNe that 
we do not use in this analysis. We compute the SNe likelihood 
using a module implemented in COSMOSIS, which re- 
produces the results in Ref. [208]. The constraint from only 
DES SNe on Qm is shown in Fig.[13|(purple). 

Galaxy clusters: The DES Y1 redMaPPer catalog consists 
of ~6500 clusters with richness larger than A = 20 in the 
redshift range z € [0.2,0.65]. The first cosmological anal- 
ysis of DES clusters (‘Clusters 1’: orange contours in 
Fig. [13). which combines cluster counts data and mass esti- 
mates from the stacked weak lensing analysis of [211], found 
a larger than 2c tension with the other DES probes. This 
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is driven by low-richness systems, and has been interpreted 
as unmodeled systematics that affect the stacked weak lens- 
ing signal of the optically selected sample. This interpreta- 
tion is supported by the analysis of (‘Clusters 2’: green 
contours in Fig. |13), which recovers results consistent with 
the other DES probes by combining cluster abundances with 
the large-scale auto-correlations of galaxy and cluster position 
and cross-correlations of cluster position with galaxy position 
and shear from DES Y1 data (4x 2pt+N). The conclusions of 
[105] are further corroborated by the analysis of which 
derive cosmological posteriors consistent with by analyz- 
ing the DES Y1 redMaPPer cluster abundances, but replacing 
the stacked weak lensing mass estimates of with multi- 
wavelength follow-up data from the SPT-SZ 2500 deg? survey 
[2.13]. 

Baryon acoustic oscillations: A sample of 7 million galax- 
ies from the DES Y3 ‘Gold’ catalog is selected in the redshift 
range 0.6 < z < 11 and used to measure the scale 
of the BAO feature in the distribution of galaxies at an effec- 
tive redshift ze = 0.835 [215]. We use a likelihood from 
Ref. for the ratio of the angular diameter distance D 4 
at Zef and the sound horizon distance at the drag epoch, ra, 
which is implemented in COSMOSIS. The simulated galaxy 
catalogs used in the analysis to derive the uncertainty of the 
measurement are described in Ref. [216]. While the BAO and 
3x2pt analyses probe common sky area and redshift range, 
and the measurements of this work include scales impacted 
by the BAO feature, the overlap in galaxy sample is small and 
the method for inferring the BAO distance ratio likelihood is 
insensitive to cosmology, so we neglect this non-zero correla- 
tion when combining the measurements. This will be further 
validated in future work that combines and studies all final 
DES Y3 probes. 

DES Year 1 3x2pt: We reanalyze the DES Y1 3x2pt data 
in the Y3 model and prior space, but do not update the scale 
cuts or marginalize over a free point mass for galaxy-galaxy 
lensing. We also make no changes in priors on systematic 
parameters (e.g., photo-z or shear calibration parameters). 

The comparison of these cosmological constraints in 
ACDM using the DES probes is shown in Fig. The com- 
bination of DES Y3 3x2pt, SNe Ia, and BAO data is also 
shown in (blue). While the constraint in Sg is driven primarily 
by 3x2pt, there is substantial gain in other parameter dimen- 
sions due to the additional data. The marginalized parameter 
values are summarized in Tables [II] & [III]and Figs. [7] & [10] 


VII. COMPARISON WITH EXTERNAL DATA 


It has been demonstrated that various combinations of low- 
redshift data and high-redshift data from the CMB can inde- 
pendently fit the ACDM model. However, the most stringent 
tests of the model will come from combining these data and 
testing whether the model can simultaneously fit the diverse 
set of cosmological probes available to us at all redshifts si- 
multaneously. These data sets are sensitive to the growth of 
density perturbations, the expansion and geometry of the Uni- 
verse, or both, and are sourced from a variety of very dif- 


ferent physical processes. The combination of the indepen- 
dent external low-redshift probes with DES Y3 data further 
reduces the potential impact of any residual systematic effects 
in the low-redshift anchor of the test, while the combined DES 
probes have been carefully calibrated from the same data and 
consistently protected against confirmation bias. Both of these 
considerations give us further confidence, for complementary 
reasons, in testing the ACDM model across the age of the Uni- 
verse. 


A. External data sets 


The likelihoods from data sets external to DES include: 

Type Ia supernovae: The Pantheon sample com- 
bines the distance measurements from 1048 SNe ranging from 
0.01 < z < 2.3, supplementing Pan-STARRS1 measure- 
ments with other available samples. 

Baryonic acoustic oscillations and redshift-space distor- 
tions: We use the constraints from SDSS measurements of 
BAO and RSD in eBOSS DR16f] When using constraints on 
f(z)og(z), where f(z) is the growth rate, from RSD, we use 
the released covariance matrices between the constraints from 
the BAO and RSD measurements. These measurements are 
expressed in terms of the Hubble distance Dy, sound hori- 
zon distance rg, comoving distance D m, and volume-average 
distance Dy = (D?,Dy,2)'/3. The measurements include 
(from low to high redshift measurements): 


* The measurement of Dy at an effective redshift of 
Zeg = 0.15 using the Main Galaxy Sample (MGS) 
[218] and adding to fog measurement from Ref. [219]. 


A re-analysed version of BOSS DR12 measurements of 
Dm and Dy from BAO, and fog from RSD, at zef = 
0.38 and 0.51 [220]. 


The eBOSS DR16 measurements of Dy /rq, Dg /rqa 
from BAO, and adding fog from the full-shape infor- 
mation, at zef = 0.698 using Luminous Red Galaxies 


(LRG) [221] ; 

The eBOSS DR16 measurements of Dy /ra when using 
BAO alone and Dy /ra4, Dy /ra, fog when using BAO 
and the full-shape information, at zep = 0.845 using 
Emission Line Galaxies (ELG) [223], 

The eBOSS DR16 measurements of Dy /ra, Dg /ra 
from BAO, and adding fog from the full-shape infor- 
mation, at Zef = 1.48 using the Quasar Sample (QSO) 
[2241 [225], 

The eBOSS DR16 measurements of Dm, Dy at zog = 
2.33 using the Lyman-o forests [226]. This data set 
only has information from BAO. 


CMB: We use the likelihoods from the Planck 2018 data 
release . Our fiducial combination of Planck likeli- 
hoods includes: 


8hhttps://svn.sdss.org/public/data/eboss/DRl6cosmo/ 
tags/v1, 0, 0/likelihoods/ 
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e The Plik likelihood of the temperature power spec- 
trum CIT in 30 < £ < 2508 and the E-mode power 
spectrum QU and the cross power-spectrum between 
temperature and E-mode C7 in the range 30 < £ < 
1996. 


* The Commander likelihood of the temperature power 
spectrum EV in 2 < £ < 29. 


e The SimA11 likelihood of the E-mode power spectrum 
CEP in2 < L< 29. 


We also use the likelihood of the lensing potential @ power 
spectrum p ? measured by Planck in the range 8 < £ < 400, 
either in combination with our fiducial combination of Planck 
likelihoods described above or alone. In the latter case we use 
the likelihood marginalized over the CMB power spectrum. 
Planck CMB will refer to the primary CMB anisotropy data 
(without lensing) unless otherwise stated. 

Big Bang nucleosynthesis: We construct an Qh? constraint 
based on observations of damped Lyman-o systems [227]. 
The primordial deuterium-to-hydrogen ratios measured from 
these systems can be translated to constraints on QyÀ? via 
Big Bang nucleosynthesis (BBN) calculations, but different 
assumptions on the BBN physics, in particular on the rate 
of the d(p, y)*He nuclear reaction, yield different final con- 
straints on Qph?. Our constraint conservatively incorporates 
the two major categories of such assumptions, namely the the- 
oretical approach presented in Ref. and the experimental 
measurement-based approach from Ref. [228]. Specifically, 
we adopt the mean and the statistical uncertainty on QA? 
from Ref. [228], and in addition introduce a systematic un- 
certainty defined by the difference between 1) the mean from 
Ref. and 2) an inverse-variance weighted average of the 
two respective means [227][228]. This results in our adopted 
constraint of 10005? = 2.195 + 0.028. 

Local Hubble parameter: We use a local h prior from 
SHOES [93], which constrains h = 0.732 + 0.013 using a lo- 
cal distance ladder that depends on measurements of Cepheids 
and type Ia supernovae. 

We use versions of these likelihoods implemented as mod- 
ules in COSMOSIS, which are used to obtain the constraints 
presented in the following. 


B. High redshift vs. low redshift in ACDM 


One of the most stringent tests of ACDM is to compare 
the prediction of the state of the Universe and amplitude of 
perturbations from the epoch of recombination, which we can 
observe from the CMB, to the current day, which we observe 
with low-redshift surveys like DES. At the time of the CMB, 
the Universe was very hot and dense, and its physics was dom- 
inated by radiation. DES is most sensitive to a period in the 
Universe approximately eight billion years later, where pertur- 
bations have grown by several orders of magnitude and non- 
linear growth is important. DES observes a volume of the Uni- 
verse spanning nearly nine billion years of its evolution. The 
volumes probed by large low-redshift surveys provide signif- 
icant additional information on potential changes to the evo- 
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FIG. 14. A comparison of marginalized constraints from three simi- 
larly constraining sets of cosmological probes in ACDM. Combined 
external BAO, RSD, and SNe Ia data (Ext. Low-z) are shown in or- 
ange, the combination of DES galaxy clustering and weak lensing 
data (3 x2pt) is shown in black, and Planck CMB (no lensing) data 
is shown in green. The three share a common parameter space in 
the Qm—Ss plane at their 68% CL bounds. The combination of Ext. 
Low-z data with DES 3x2pt is shown in purple and this combined 
additionally with Planck CMB (w/ lensing) is shown in blue. 


lution of perturbations or growth of the Universe over time, 
allowing them to strongly test the nature of dark energy. 

By taking precise measurements of the ACDM model from 
CMB observations and predicting what we should observe in 
terms of the amplitude of perturbations or matter density in 
the late Universe, we can test whether our observations from 
surveys like DES agree with those predictions. If they do not 
agree at high significance, we have demonstrated that ACDM 
cannot describe the full evolution of the Universe. There has 
been considerable debate about the tendency of late Universe 
measurements to prefer slightly lower matter density or am- 
plitude of clustering relative to measurements from the CMB 
(e.g., [131271 B3][101] [102]). As more powerful data becomes 
available, like the current DES Y3 analysis, we can determine 
whether these measurements converge towards or away from 
the Planck CMB prediction. 
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FIG. 15. A comparison of marginalized parameter constraints from 
three similarly constraining sets of cosmological probes in wCDM. 
Combined external BAO, RSD, and SNe Ia data (Ext. Low-z) are 
shown in orange, the combination of DES galaxy clustering and 
weak lensing data (3x2pt) is shown in black, and Planck CMB (no 
lensing) data is shown in green. The combination of Ext. Low-z data 
with DES 3x2pt is shown in purple and this combined additionally 
with Planck CMB (w/ lensing) is shown in blue. 


We compare three similarly constraining and complemen- 
tary subsets of available cosmological probes in ACDM and 
wCDM in Figs. [14] & [15] The external low-redshift SNe Ia, 
BAO, and RSD data primarily constrain Qm and w, while 
the DES 3x2pt data adds substantial information on A, or 
og, Which helps to further constrain Qm and w through 
degeneracy-breaking of correlated parameters. These external 
low-redshift data sets complement the DES weak lensing and 
large-scale structure information by probing the growth and 
geometry of the cosmological model in fundamentally differ- 
ent ways. The CMB is able to tightly constrain both Qm and 
Ag or og in ACDM, but is comparable in constraining power 
to DES 3x2pt in wCDM, since it primarily has access to in- 
formation limited to the surface of last scattering at z ~ 1100. 
The combination of DES 3x2pt with the other low-redshift 
data provides substantial gain in As, og, Qm, and w. We list 
marginalized parameter constraints for these probes in Tables 
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FIG. 16. A comparison of the marginalized parameter constraints 
in the ACDM model from the Dark Energy Survey with predictions 
from Planck CMB data (no lensing; green). We show the fiducial 
3x 2pt (solid black) and the combined Y3 3 x 2pt and Planck (orange) 
results. 
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1. Consistency results 


We show the comparison of DES 3x2pt and the Planck 
CMB data for the ACDM and wCDM models in Figs. 
& Visually, we find better agreement in the overlap of 
the marginalized Qm—Sg parameters with the DES Y3 3x2pt 
data than found in the DES Y1 analysis [33], despite substan- 
tial improvements to the precision of both DES and Planck 
predictions. This is qualitatively unchanged when using the 
more precise, optimized ACDM version of the analysis that 
uses more small scale information — the DES contour shrinks, 
but asymmetrically in the direction of the CMB prediction. 

We evaluate the consistency of the DES and Planck data 
in several ways, including shifts in parameter space and the 
Bayesian evidence. These are described further in Sec. [IV E] 
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FIG. 17. A comparison of the marginalized parameter constraints 


in the wCDM model from the Dark Energy Survey with predictions 
from Planck CMB data (no lensing; green). We show the fiducial 
3x 2pt (solid black) and the combined Y3 3 x 2pt and Planck (orange) 
results. 


and full results are provided in App. |F| We find a parame- 
ter difference of 1.50 (p = 0.13) in the cosmological model 
space and a Suspiciousness of 0.7 + 0.1c, corresponding to 
p = 0.48 + 0.08. This generally leads to the conclusion 
that despite substantially increased precision from both ex- 
periments, we find no significant evidence against the ACDM 
model from comparing these data sets. Agreement between 
DES and Planck in these metrics has improved relative to 
the comparison of DES Y1 3x2pt and earlier Planck results, 
which gave a parameter difference of 2.2c and Suspiciousness 
of 2.4 + 0.20 [133]. The combined DES and Planck CMB 
contour is shown in orange in Figs. [16] & 


We repeat this exercise for the full combined low-redshift 
data, including DES 3 x 2pt, all BAO, and external SNe Ia and 
RSD data. This comparison is shown in Figs. [14] & [15] and 
is highly complementary, as the external probes are sensitive 
to both growth and geometry in the model in ways the DES 
3x2pt data is not, and come from a variety of different exper- 


iments. We find better agreement between all of these low- 
redshift probes and Planck CMB predictions than in the com- 
parison with DES 3 x2pt data alone, with a parameter differ- 
ence of 0.90 or p = 0.34. These results indicate that we can 
combine all these available cosmic probes into a single joint 
result in the following subsection. 


There are several reasonable motivations for caution in the 
interpretation of any strong evidence for or against cosmolog- 
ical consistency in tests like this. It is worth noting that while 
we have multiple redundant low-redshift sources of informa- 
tion for each main cosmological probe used, it would be useful 
to have a second, blinded large-scale CMB polarization ex- 
periment to increase confidence in the test at the high-z limit. 
While polarization data is required to break degeneracies in 
the cosmological parameters with the optical depth 7, we also 
repeat the caution from Ref. against over-interpreting the 
Planck polarization results and the sensitivity of the final pa- 
rameter constraints to assumptions made in the construction 
of the likelihood, which can lead to a «1o shift toward the 
DES posterior relative to the fiducial Planck likelihood. Sim- 
ilar shifts are seen based on certain analysis choices in the 
DES results as well, which are shown in App. |E| Neither the 
shift in Planck posteriors or those from other analysis choices 
in DES contribute to a significant change in the final inter- 
pretation of the comparisons. Finally, the DES Y3 analysis 
has uncovered potential systematics connected to photometry 
(e.g., Sec.|V C). While there is evidence that these do not im- 
pact the cosmological results, and thus would not impact this 
comparison of data sets, they have not been connected to a 
specific source. However, these are examples of unresolved 
uncertainties that call for additional care in interpreting any 
statements about the consistency of early- and late-universe 
probes in ACDM, which should be addressed for future more 
precise analyses. 


C. Joint cosmological constraints in ACDM and wCDM 


We find that external low-redshift (BAO+RSD+SNe Ia), 
Planck CMB, and DES 3x2pt data sets are able to provide 
three independent, highly complementary, and similarly pow- 
erful constraints on parameters related to dark matter and dark 
energy in the ACDM and wCDM models, as seen in Figs. 
& Given the results of the above consistency tests, de- 
tailed in App. |F| these data sets are each consistent with one 
another, and thus can be combined into a joint constraint on 
the models. We present these joint results in Figs.[14]&[15]and 
a summary in Figs. [7] & [TO|and Tables [IT] & [III] In the ACDM 
model, we find 


Sg —0.812*0 095 (0.815) 


Om = 0.3067) 005. (0:306) Q2) 


og = 0.804* 0-005. (0.807). 
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In the wCDM model, 


og = 0.810*0 010. (0.804), 
Nm 50.309 79106. (0.208), (23) 
w = —1.031*0029 (—1.001). 


We find R = 7.8, indicating that there is also no preference 
for wCDM over ACDM in the full joint data analysis. 


These data sets together are able to provide unprecedented 
precision on the cosmological parameters of the models. In 
ACDM, we are able to constrain og, Sg, h, Qp, and n, to less 
than 196; Qm and A, to about 1%; 7 to about 10%; and place 
an upper limit on the sum of neutrino masses of 5 m, < 0.13 
eV (95% CL). In wCDM, we are able to constrain n, to less 
than 196; Qm, Qp, h, and A, to about 1-296; w to about 3%; T 
to about 10%; and place an upper limit on the sum of neutrino 
masses of $5 m, < 0.17 eV (95% CL). Individually, the three 
subsets of data constrain og and Qm in ACDM with FoM be- 
tween 2000 and 4000, while combined, they reach a FoM of 
34,000. This clearly demonstrates the highly complementary 
nature of these three independent data sets. 


D. Comparison of Lensing Probes 


We are able to probe the distribution of large-scale structure 
via weak gravitational lensing (cosmic shear) in two very dif- 
ferent ways, using either the shapes of galaxies or the CMB 
photons as tracers for the reconstruction of deflections in the 
path of the light. These probe the same physical phenomenon 
via independent sources and measurement methods, which are 
sensitive to different types of systematics. While the effective 
kernel of CMB lensing is sensitive to higher redshift 
structure than galaxy lensing, their comparison provides a sig- 
nificant validation of the robustness of modern weak lensing 
results. We show this comparison in Fig. where we find 
very good agreement between the two cosmic shear measure- 
ments and the full 3x2pt measurement from DES. 

In addition to DES, other concurrent photometric surveys 
HSC and KiDS are also pursuing preci- 
sion weak lensing measurements using galaxy shapes. These 
three surveys span a range of depth and survey area tradeoffs, 
with HSC being deepest, DES widest, and KiDS using the 
widest wavelength coverage. We over-plot recent results from 
each of the surveys with our DES cosmic shear and 3x2pt 
results in Fig. Unlike other comparisons in this work, 
these external survey data have not been re-analyzed within 
a consistent model and prior space. Thus, no direct or rig- 
orous comparison can be made about data consistency. We 
defer a detailed discussion of the consistency of concurrent 
photometric weak lensing surveys and their combination (e.g., 
Ref. (230]) to a future work. The apparent orthogonal direc- 
tion of the KiDS+BOSS+2dFLenS 3 x 2pt contours to the DES 
3x2pt contours is driven by the very strong constraint coming 
from spectroscopic clustering, similar to the orientation of the 
DES y++w(0) constraint. 
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FIG. 18. A comparison of weak lensing constraints on the ACDM 
model. Weak lensing of the CMB is shown in green, weak lensing 
of galaxies in DES is shown in blue, and the combined DES 3x2pt 
data is shown in black. 
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FIG. 19. The DES ACDM-optimized 3 x 2pt and cosmic shear, HSC 
and KiDS cosmic shear, and KiDS lensing + BOSS+2dFLenS spec- 
troscopic 3x 2pt data results are over-plotted for the ACDM model. 
Unlike other comparisons in this work, these external survey data 
have not been re-analyzed within a consistent model and prior space. 
Thus, no direct or rigorous comparison can be made about data con- 
sistency. 


26 


| [ACDM — Fid. 3x2pt+BAO+BBN 
0.38 F == DES 3x2pt+All Ext. 
——  BAO--BBN 
, — Planck CMB 
Msn = SHOES Local h 
0.34L 
€ 932r 
0.30} 
0.281 
0.26. 


0.64 0.66 0.68 0.70 0.72 0.74 0.76 
h 


FIG. 20. Marginalized constraints on h and Qm in the ACDM model 
are compared to the SHOES local determination of h. Planck CMB 
data and the combination of BAO and BBN data provide compara- 
ble uncertainties on h compared to the local constraint. Adding DES 
3x2pt to BAO and BBN improves the constraint on A slightly due 
to 3x2pt providing additional information on Qm, while the combi- 
nation of DES 3 x2pt and all non-local external data provide a con- 
straint on A that is a factor of 3-4 more powerful than the local deter- 
mination. 


E. Constraints on the Hubble parameter 


There is an interesting disagreement in local measurement 
of the Hubble parameter h and marginalized constraints on 
h from cosmological experiments. Multiple local measure- 
ments prefer a higher value of the expansion velocity, such 
as, most prominently, the astronomical distance ladder (e.g., 
h = 0.732 + 0.013 with Cepheid variable stars; h = 
0.733 + 0.040 with Mira variable stars), or masers (e.g., 
h = 0.739 + 0.030km/s/Mpc [95]). These local mea- 
surements stand in contrast to constraints from the CMB by 
Planck, which prefer h = 0.6655 Dus (when the neutrino 
mass density is varied) [13]. However, there are also local 
measurements with lower values reported (h = 0.696 + 0.019 
with tip of the red giant branch distance ladder; h — 
0.074 * 505. with strong lensing when combining the TD- 
COSMO+SLACS data set [97]). The Hubble tension may 
indicate new physics and it is crucial to improve measure- 
ments, revisit assumptions [e.g., [98], check for consis- 
tencies among different measurements, and invest in novel, 
independent methods and probes [99] [100]. 


We can also constrain the value of h independently of CMB 
data using a combination of BAO, BBN constraints on Qyh?, 
and DES 3x2pt measurements. Constraints on h and Qm in 
ACDM are summarized in Fig. The determination of h 
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FIG. 21. Marginalized constraints on the sum of neutrino masses 
in the ACDM model. We show the DES fiducial 3 x 2pt constraints 
(black), DES 3x 2pt combined with external BAO, RSD, and SNe 
Ia (orange), Planck CMB constraints (green), and DES 3 x 2pt com- 
bined with all of these external data sets. The upper panel shows 
the one-dimensional marginalized posteriors for X` mz, with shaded 
95% confidence regions. The lower panel shows 68 and 95% CL for 
Om and 35 mp. 


using BAO and BBN is of similar constraining power to that 
of the CMB and agrees very well with the CMB constraint on 
h. Adding DES 3x2pt data slightly improves the constraint 
on h and shifts it to higher values by about 1c. Combining 
DES 3x 2pt data with BAO, RSD, SNe Ia, and Planck CMB 
(w/ lensing) leads to a marginalized constraint on h 


h = 0.680*0005 (0.681) (24) 


that is 3-4 times more powerful than any current local mea- 
surement of h. Constraints on other cosmological parameters 
are summarized in Tables[II]& [M] We find no significant im- 
pact on the other cosmological parameters by adopting this 
high-redshift anchor for the expansion rate vs a local prior on 
the expansion rate from Ref. [93]. The final joint constraint 
on h is consistent with the Planck- or BAO+BBN-only con- 
straints and slightly less than 4c offset relative to the local h 
by SHOES. 


F. Neutrino Mass 


Figure shows marginalized constraints on the sum of 
neutrino masses, where neutrino mass and density €), are re- 
lated via 35 m, = 93.140, h? eV. We model massive neutri- 
nos as three degenerate species of equal mass. As expected, 
DES does not constrain neutrino mass: whether alone or in 
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combination with external BAO, RSD, and supernova data, 
the marginalized posterior of the sum of neutrino masses is 
bounded by its prior. Given this, DES data is expected to add 
very little direct information on neutrinos. 

Beyond constraints on neutrino mass themselves, a motiva- 
tion for looking at the neutrino mass constraint is to highlight 
a feature of the relationship between Planck and low-z con- 
straints on Qm, as shown in e.g. Fig. [14] As has been pre- 
viously discussed [183], a geometric degeneracy means that 
CMB-only constraints are unable to distinguish between Qm 
and 5 m,, but combining CMB data with low redshift expan- 
sion history constraints from BAO and/or SNe can break that 
degeneracy. This also illustrates that the parts of the Planck 
CMB posterior at higher Qm also have relatively high neutrino 
mass. This is further supported by the behavior of constraints 
when neutrino mass is fixed, as shown in App.[E] 

Lowering the clustering amplitude has a similar impact on 
Qm, which we see as an increase in the upper limit on $5 m, 
of 23% when combining DES 3 x2pt with Planck. Combin- 
ing the DES 3x2pt data, other low-redshift data, and Planck 
CMB, we find an upper limit 


Y m, < 0.13 eV (95% CL), (25) 


nearly a factor of three reduction from the CMB-only con- 
straint. 


VIII. CONCLUSIONS 


We have described the 3x2pt measurements, calibration, 
modeling, and analysis from the first three years of DES data. 
The substantial improvement in statistical power of the DES 
Y3 data, which cover an area of sky about three times that of 
DES Y1, has required substantial improvements in almost ev- 
ery part of the data processing, analysis and inference. The 
specific improvements relative to the DES Y1 analysis are 
detailed in App. |B| but we also briefly summarize here. At 
the catalog level, some of the most important updates for this 
analysis include improved PSF modeling and two comple- 
mentary lens sample selections. We have substantially revised 
our shear and redshift inference and calibration processes. 
This includes more realistic image simulations to derive cor- 
rections on shear and redshift bias due to blending and detec- 
tion, and a redshift inference process that combines spectro- 
scopic and deep, multi-band photometric redshifts from DES 
deep-field data, cross-clustering between source and high- 
quality photo-z and spectroscopic samples, and small-scale 
galaxy-galaxy lensing shear ratio information. Finally, we 
have updated several components of the analysis, including 
blinding, modeling of non-Limber and RSD contributions to 
the clustering signal, mitigation of nonlocal effects in y+, im- 
proved validation of the covariance matrix, and improved, cal- 
ibrated metrics for evaluating the internal and external consis- 
tency of data sets. 

The statistical power of the DES Year 3 data set has posed 
unique challenges for precision cosmological inference, some 
of which were unforeseen. We have identified some puzzling 


results from our photometric lens samples and galaxy cluster- 
ing, which have not been identified with a clear source. This 
can be seen at very large significance as an apparent disagree- 
ment in the clustering and lensing amplitudes at all redshifts 
and angular scales for redMaGiC, and in the highest redshifts 
of MagLim not used in the fiducial analysis. After unblinding 
our analysis, it was necessary to make two important revisions 
to the fiducial analysis plan: 1) we made the MagLim galaxy 
sample our fiducial lens sample owing to the decorrelation of 
the lensing and clustering amplitudes inferred from the red- 
MaGiC sample, and 2) we dropped the two highest redshift 
bins of the MagLim sample as they contributed to a very poor 
fit to all models considered in this paper. Sections [VA]& 
provide the detailed rationale. These decisions were made af- 
ter extensive, careful investigations of possible systematics, 
which did not reveal problems that we could address. Further 
investigations are already underway and may reveal the source 
of these issues, but our robustness tests so far indicate that any 
potential changes to the results will lie well within our quoted 
uncertainties. 

We have achieved a factor of two improvement in statis- 
tical power relative to the DES Y1 3x2pt analysis in the 
og—Qm marginalized parameter plane, providing competitive 
cosmological constraints relative to both the combination of 
all other external non-lensing low-redshift data and the pre- 
dictions from the Planck CMB data. We find consistent 
cosmological constraints from cosmic shear and the com- 
bination of galaxy clustering and galaxy—galaxy lensing, as 
well as consistent cosmological results from 3x 2pt utilizing 
two lens samples, MagLim and redMaGiC. For the fiducial 
3x 2pt analysis in ACDM, we find constraints on the cluster- 
ing amplitude Sg = 0.776+9:017 (0.776) and matter density 
Qn = 0.339+9:93? (0.372). In the wCDM model, we find 
Om = 0.352+9:933 (0.339), and dark energy equation of state 
parameter w = —0.98*0:22 (—1.03). 

The low-redshift measurements of the matter clustering am- 
plitude by some galaxy surveys have tended to find lower vari- 
ance relative to the prediction from Planck CMB anisotropies, 
which may indicate some inconsistency between the low- and 
high-redshift Universe within the ACDM model, with claims 
of up to 2-30 significance. The DES Y3 3x2pt analysis 
is an ideal experiment to test whether this is a real problem 
with the ACDM model. There have been substantial improve- 
ments in constraining power for both DES and Planck since 
the DES Y1 3x2pt analysis, yet we continue to find that DES 
3x2pt data and the combination of 3x2pt with BAO and ex- 
ternal SNe Ia and RSD data that the low-redshift Universe are 
consistent with predictions in the ACDM model from mea- 
surements at the time of the CMB from Planck. We find all 
three independent data set combinations (DES 3x2pt; BAO, 
RSD, and SNe Ia; and Planck CMB) to be mutually consistent 
within ACDM. Despite caveats on the precision with which 
we can make this statement discussed in Sec. this is the 
most powerful test of the standard cosmological model to date, 
comparing predictions from measurements of acoustic peaks 
in the early plasma of the Universe when it was 380,000 years 
old to measurements of large-scale structure from low-redshift 
surveys like DES spanning nearly nine billion years of cosmic 
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evolution to the current day. 

Combining DES 3x2pt, CMB, BAO, RSD, and SNe 
Ia data allows us to place the most precise con- 
straints on the ACDM and wCDM models to date. 
We find Sg = 0.81273 008 (0.815) and Qn = 
0.306*0002 (0.306) in ACDM; og = 0.812* 0005 (0.804), 
Qm = 0.302*0005 (0.298), and w = —1.031*0 03» (—1.00) 
in wCDM. Additionally, we find an independent constraint 
on the Hubble parameter combining DES 3x2pt, BAO, and 
BBN data of h = 0.676*0 009 (0.673), which is consistent 
with the Planck prediction for h. From the combination of 
DES 3x2pt, CMB, BAO, RSD, and SNe Ia data, we find 
h = 0.680*0004 (0.681). This is slightly closer to the lo- 
cal h measurement by SHOES than Planck, but a factor of 
three to four more constraining than either the local or Planck 
measurement of h. We are also able to constrain the sum of 
neutrino masses to be $5 m, < 0.13 eV (95% CL) in ACDM. 

While we have shown that the inferred cosmological con- 
straints from the fiducial analysis are robust, there remains 
significant work to fully characterize the underlying causes of 
these potential systematics and examine other potential theo- 
retical causes in extended model spaces beyond the dark en- 
ergy models considered in this work. Further understanding 
the potential systematic issues related to differences between 
photometric clustering and galaxy-galaxy lensing; improve- 
ments in how we deal with shear calibration and redshift infer- 
ence in the presence of blending; and finding ways to improve 
systematic floors in our redshift inference are all important 
next steps. This continued followup work will be critical to 
the final DES Year 6 analyses and future “Stage IV’ photo- 
metric surveys like the Euclid Space Telescope] the Nancy 
G. Roman Space Telescope] "|and the Vera C. Rubin Obser- 
vatory Legacy Survey of Space and Time (LssT)["] 

The novel advances required for the DES Y3 analyses, sum- 
marized in the accompanying 29 papers [11511143], set the 
stage for these future precision low-redshift large-scale struc- 
ture and weak lensing studies. DES has utilized only half its 
final data set in the DES Y3 3x2pt and BAO analyses, and 
future SNe Ia and galaxy cluster analyses promise even larger 
improvements in statistical power. The legacy analyses of the 
full DES data will be a focus of the next several years leading 
up to the start of Stage IV dark energy experiments. 
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Appendix A: Summary of Associated Papers 


This paper is built on the results presented in 29 accompa- 
nying papers. A useful way to navigate these is to divide them 
up into five categories: 

Catalog Papers: The link between the raw images and the 
two-point functions from which cosmology is extracted is a 
set of catalogs. The Gold catalog (Sevilla-Noarbe et al. [115]) 
uses coadd images to identify galaxies and their properties. 
This is a first step to almost all ensuing work. Estimating the 
redshifts of those galaxies hinges in large part on the much 
deeper catalog of galaxies from the DES deep fields with over- 
lapping near-infrared photometry presented in Hartley, Choi 
et al. [116]. Finally, the shear catalog is presented in Gatti, 
Sheldon et al. [117], which uses a new PSF measurement de- 
scribed in Jarvis et al. [118]. 

Simulations: To test our models and calibration of the data, 
we rely on large suites of cosmological and image-level sim- 
ulations. The Buzzard simulations (DeRose et al. [119]) are 
generated from cosmological N-Body simulations, populated 
with realistic galaxy samples, and then used in end-to-end 
analyses to stress test our modeling and methods. To cal- 
ibrate blending and detection biases, we generated multiple 
image simulations, described in MacCrann et al. [120], where 
shear and redshift biases are evaluated. Balrog (Everett et al. 
([121]) is a tool developed to inject realistic images of galaxies 
into real DES images to evaluate the survey selection function. 
This is critical for multiple purposes, including our photomet- 
ric redshift inference. The full DES detection and measure- 
ment pipelines are run on the simulations described in Mac- 
Crann et al. and Everett et al. [121]. 

Photometric Redshifts: The redshift distributions for the 
redMaGiC and MagLim lens galaxy samples are validated us- 
ing cross-correlations with spectroscopic galaxies in Cawthon 
et al. [122], which informs our priors on the uncertainty of 
the redshift distribution. The redshift inference process for 
the source galaxies is a much more involved process. The 
overview of this work is presented in Myles, Alarcon et al. 
[123], which builds on the self-organizing map formalism 
developed in Buchs, Davis et al. [124], and incorporates 
constraints from the cross-correlation with both redMaGiC 
and spectroscopic galaxies described in Gatti, Giannini et al. 
[125]. A followup analysis applying this methodology to the 
MagLim lens redshift distributions is described in Giannini 
et al. [126]. The sampling of the resulting realizations of the 
redshift distribution are described in Cordero, Harrison et al. 
[127]. Finally, we also add additional information from the 
posterior from the shear ratio likelihood, which is most sen- 
sitive to photo-z parameters and is described in Sánchez, Prat 

Analysis: The modeling outlined in the text is described and 
tested in detail in Krause et al. [129]. One of the most impor- 
tant pieces of the analysis is the generation of the covariance 
matrix of the two-point functions. The way in which this is 
generated and tested is described in Friedrich et al. [130]. As 
mentioned in the text, we remain blinded to the data through- 
out the analysis to avoid unconscious bias; this blinding has 
multiple levels, but one of the key tiers is described in Muir 


et al. [131]. Methods to evaluate the internal consistency of 
our data and to evaluate consistency with external data sets 
were calibrated and described in Doux et al. & Lemos 
et al. [133]. Efficient sampling of the likelihood can save 
months of analysis time, and ways to optimally, and more im- 
portantly robustly, sample our likelihood are summarized in 
Lemos et al. (134]. Finally, the optimization of the MagLim 
lens sample is described in Porredon et al. [135]. 

Results: The galaxy clustering weights and measurement 
are presented in Rodríguez-Monroy et al. [136]; galaxy- 
galaxy lensing in Prat et al. and in Sánchez, Prat et al. 
for the smaller-scale shear ratio measurements; and the 
weak lensing convergence mass map in Jeffrey, Gatti et al. 
[138]. The results of combining clustering with galaxy- 
galaxy lensing are presented in three papers, one of which fo- 
cuses on magnification (Elvin-Poole, MacCrann et al. [139]), 
another on redMaGiC and the bias model (Pandey et al. 
[140]), and the last on the MagLim sample (Porredon et al. 
[141]. Finally, the cosmic shear results are presented in two 
papers, one of which focuses on observational systematics 
(Amon et al. (142]) and the other on biases in theoretical mod- 
eling (Secco, Samuroff et al. (143]). 


Appendix B: Differences relative to DES Year 1 analysis 


In this Appendix, we summarize the major differences in 
the DES Y3 analysis relative to DES Y1. 

Data Processing: DES Y3 contains significantly more data 
than DES Y1 [144], which are processed with an improved 
version of the DESDM system [146]. Photometric calibra- 
tion is performed with the forward global calibration method 
[45], which significantly improves the relative photometric 
calibration compared to the calibration techniques applied in 
DES Y1 [233]. The Y3 catalogs also introduce more robust 
morphological classification based on multi-epoch fitting and 
improved flagging to enhance the quality of the galaxy sam- 
ple. Further details on improvements to the DES data set and 
object catalogs can be found in Sevilla-Noarbe et al. (115). 

Catalog-level: We have produced a calibrated deep-field 
data reduction for use in the analysis. We have constructed 
only one shape catalog (METACALIBRATION; [117]), but pro- 
duced two very different lens samples to compare cosmolog- 
ical results [136] [140] [141]. We use a new PSF model (PIFF) 
[118]. 

Simulation & Calibration: We have produced the Bal- 
rog simulation to characterize the wide-field survey selection 
function of objects derived from the deep-field coadd images 
[21]. We have produced a new suite of image simulations for 
shear calibration that are more realistic and better match the 
data [120]. We have produced a new methodology for shear 
calibration that explicitly accounts for the impact of blend- 
ing and detection biases as a function of redshift, including 
modifications to the effective n(z) [120]. We have developed 
a new redshift inference and calibration framework, mixing 
spectroscopic and photometric redshift information, to pro- 
duce Bayesian posterior n(z) samples [12211128]. 

Modeling: We have improved the rigor of constraining po- 
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TABLE IV. The parameter differences for the redMaGiC analysis 
relative to Table[I] 


Lens Galaxy Bias (redMaGiC ) 


bili € [1,5]) Flat (0.8, 3.0) 
Lens magnification (redMaGiC ) 

Ci Fixed 0.63 

c? Fixed —3.04 

C? Fixed —1.33 

ci Fixed 2.50 

CP Fixed 1.93 
Lens photo-z (redMaGiC ) 

An Gaussian (0.006, 0.004) 
Az? Gaussian (0.001, 0.003) 
Az Gaussian (0.004, 0.003) 
Az} Gaussian (—0.002, 0.005) 
Az Gaussian (—0.007, 0.010) 
ot (i € [1,4]) Fixed 1.0 

021 Gaussian (1.23, 0.054) 


tential bias from model approximations, with explicit accu- 
racy goals at the x? and parameter levels [129]. We have in- 
cluded non-Limber and redshift-space distortion (RSD) con- 
tributions to the clustering theory [129]. We have included 
the impact of lens magnification [139]. We have made 
updates to the covariance modeling and improved the rigor 
of validation tests in both x? and parameter space [130]. We 
utilize a new intrinsic alignment model that allows for 'red', 
*blue', and mixed alignment modes to second order in pertur- 
bation theory [174]. We have developed a nonlinear galaxy 
bias model for the analysis [140]. We have updated our 
methodology for determining scale cuts using explicit accu- 
racy goals in both x? and parameter space [129]. 

Analysis: We have utilized a new summary-statistic-level 
blinding scheme [131]. We have introduced a new shear- 
ratio likelihood [128]. We have introduced a new cross- 
clustering redshift likelihood [125]. We have imple- 
mented a method for sampling over the full-shape n(z) sam- 
ples [127]. We analytically marginalize over uncertainty in 
the lens sample clustering weights [136]. We analytically 
marginalize over a point-mass contribution to ^, to mitigate 
non-local effects [129]. We marginalize over the widths of the 
lens n(z) (122) [136]. We have updated and calibrated a new 
set of internal and external consistency tests [133]. We 
have updated requirements on sampling and evidence preci- 
sion and now use the PolyChord sampler [134]. 


Appendix C: Differences in redMaGiC analysis 


The redMaGiC lens analysis differs in several ways due to 
the different redshift range, binning, and sample selection rel- 
ative to the MagLim lens analysis. The parameterization and 
prior changes relative to the fiducial analysis are listed in Ta- 
ble[IV] In particular, we do not vary the width of the redshift 
distributions in the first four redMaGiC lens bins, since con- 
straints from clustering on the n(z) agree well with the predic- 
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FIG. 22. The measured w(0) correlation functions for each tomographic bin combination used in the redMaGiC analysis, which is indicated 
by the 2, j label in each set of panels. The best-fit ACDM model from the analysis using redMaGiC is plotted as the solid line in the top 
part of each panel, while the bottom parts of each panel shows the fractional difference between the measurements and the model prediction, 
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50). The best-fit model with fixed Xiens is shown in black, while the best-fit model marginalizing over 


Xiens is shown in blue. Both the top and bottom part of each panel includes 1c error bars. Small angular scales where the linear galaxy bias 
assumption breaking down are not used in the cosmological analysis; these scales are indicated by grey shading. 
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FIG. 23. The measured y (0) correlation functions for each tomographic bin combination using the redMaGiC sample, with labels as described 
in Fig. B] The best-fit ACDM model from the analysis with fixed Xie; is plotted as the solid line in the top part of each panel, with dotted 
curves indicating a negative model fit. The best-fit model marginalizing over Xjens is shown in blue. The bottom part of each panel shows 


the fractional differences between the measurements and the model prediction, (y? 


obs. 


u a") ony (with y-axis range +50). In both panels, 


lo error bars are included. Angular scales not used in the cosmological analysis are indicated by grey shading, which are excluded on small 


scales where the linear galaxy bias assumption breaks down. 
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FIG. 24. A test of the convergence of the theoretical 3 x 2pt covari- 
ance, showing marginalized parameter constraints using two itera- 
tions of the covariance that use the best-fit cosmological parameters 
from the previous analysis. 


tions from the redMaGiC photo-z algorithm [122]. The mag- 
nification parameters are also fixed to different values based 
on measurements from Balrog [139]. We also show the red- 
MaGiC galaxy clustering and galaxy-galaxy lensing data vec- 
tors and best-fit cosmological model in Figs. [22]& 3] 


Appendix D: Observer Bias & Validation Process 


The role of the observer in determining how to proceed 
in an analysis or measurement and when to accept a result 
as final has been demonstrated to contribute to uncharacter- 
ized bias in results (e.g., [235]. To protect against 
this, all cosmologically relevant measurements, calibrations, 
model validation, and fiducial analysis plans were performed 
blinded. This was done in three stages, to provide redundant 
protection: 


1. The ellipticity in the shape catalog was blinded by a 
random factor in b € [0.9, 1.1], transforming ni — bni, 
where 7; = 2e; arctanh(e)/e and e? = e? + e2. The 
factor b is generated in a random way by hashing a short 
string phrase, and this phrase was known only to the 
people making the catalog. 


2. The two-point correlation functions were coherently 
shifted by an unknown vector in cosmological parame- 
ter space (Qm, w), as described in Ref. [131]. This pro- 
duces a new data vector corresponding to an unknown 
wCDM cosmology, and has the benefit of leaving all 
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three two-point correlation functions internally consis- 
tent. 


3. All parameter values and the axes of posterior plots in- 
ferred from chains were randomly shifted. 


These protections were removed one by one as the analysis 
matured. 

First, the catalog-level blinding was removed when the val- 
idation and planned calibration of the catalog were finalized, 
to allow null tests and non-cosmological measurements to be 
remade for the papers while the data vectors were still blinded. 
In practice, this also only occurred after model validation tests 
had been completed and other analysis and calibration plans 
were also finalized. 

Once the data vector, shear and redshift calibration, all 
model validation and choices, and plans for testing internal 
and external consistency were finalized, we proceeded with a 
pre-determined set of tests on the blinded results. The enumer- 
ated tests below were those that would have led us to recon- 
sider whether to proceed with unblinding and pursue potential 
systematic causes. In addition to this, other validation checks 
were also performed at this point on the blinded data vector 
and posteriors, but were not required for unblinding. 


1. Verify the final independent redshift posteriors are con- 
sistent (i.e., from SOMPZ, the shear ratio, and cluster- 
ing cross-correlations). 


2. Verify no posteriors of systematic parameters concen- 
trate at the edge of their priors in ways that are not un- 
derstood. 


3. Verify the goodness-of-fit for each subset of the 
data (i.e., cosmic shear and galaxy clustering+galaxy— 
galaxy lensing). We used the posterior predictive distri- 
bution (PPD) with a quantitative requirement p > 0.01. 


4. Verify that cosmic shear and galaxy clustering+galaxy— 
galaxy lensing are consistent with each other, with the 
same quantitative requirement. 


5. If any of the previous tests failed in ACDM, a passing 
condition in wCDM would also be sufficient. 


The theoretical covariance matrix associated with the 
3x2pt data is calculated with an assumption about the true 
cosmological and other parameters in our model that may be 
different from what we find after unblinding our data. Using 
a covariance matrix that assumes the wrong cosmology can 
bias the cosmological inference process. To mitigate this, we 
recalculate the covariance matrix at the best-fit cosmology of 
the initial 3 x 2pt analysis and run all final chains with this co- 
variance. We confirm that this process has converged by com- 
paring the result of this analysis with an analysis based on yet 
a third covariance calculated at the second best-fit cosmology. 
This is shown in Fig. 4] 

We unblinded the redMaGiC sample before the MagLim 
sample, finding it to pass all unblinding checks. After un- 
blinding it became apparent from additional internal consis- 
tency (PPD) tests that we had statistically significant evidence 


for a potential systematic error in the clustering part of the 
data vector. This was explored at length in parallel to the final 
validation of the MagLim sample leading up to its unblind- 
ing. We verified, though did not require, that there was no 
evidence that MagLim suffered from the same effect before 
unblinding it. The MagLim sample, however, did fail the X? 
criterion above for both models after updating the parameter 
values for the covariance matrix after the initial unblinding 
chains, with an excess x? œ% 100. This was not seen before 
unblinding due partly to a poor choice of galaxy bias values 
for the initial covariance, leading to the updated covariance 
matrix being 40-5096 smaller in the w(0) block. 


We had agreed that it was implausible that our data should 
be able to distinguish a non-wCDM model at such a large x? 
before unblinding, and so we had planned to explore poten- 
tial causes for the poor model fit before proceeding in such a 
case, correcting any problems that were then discovered be- 
fore unblinding. Unfortunately, this effect was hidden due to 
the initial covariance before parameter updates being weaker 
in the w(0) blocks, and so this was a post-unblinding cor- 
rection. However, the solution (imposing an upper redshift 
cut for the MagLim sample) is something we would almost 
surely have pursued before unblinding had this problem been 
apparent then. These lens sample issues and related tests are 


summarized further in Sec.[V C]and in Refs. [136 : 


The list of changes made to the analyses post-unblinding 
are as follows. For redMaGiC, we investigated how the lens 
weights were calculated due to indications that an alterna- 
tive method that used a principal component (PC) basis of 
the observing condition maps produced a non-trivial shift in 
the clustering signal in the direction to correct the observed 
excess clustering parameterized by Xjens. However, it was 
determined that this basis was contaminated by true large- 
scale structure modes. A simplified basis limited to the first 50 
PCs to remove any contaminated modes was used in the final 
analysis that gives consistent results with the original weights 
used at unblinding. An additive component was added to 
the w(0) covariance block of both lens samples that accounts 
for potential over-correction and differences between the two 
weights methods. The MagLim sample used this final weight- 
ing and covariance when unblinded. As discussed previously, 
the highest two redshift bins of the MagLim sample were also 
removed to resolve a very poor x? fit to any dark energy model 
considered and a strong indication of inconsistency between 
w(0) and ^; approaching Xjens = 0.55 in the highest red- 
shift bin. There was no indication of the impact of Xjecns 
for MagLim generally within the redshift range overlapping 
redMaGiC. This change in redshift limit resulted in a shift 
of 0.780 in the Qm—Ss plane, almost fully in the direction 
of increasing Qm, relative to the parameter values found at 
unblinding. This is consistent with shifts observed when re- 
moving parts of the redshift range of the data, so not clearly 
evidence of systematic impact on the parameter values. 


37 


— Fid. 3x2pt 

0.84f ACUM === Hyperrank 
-— No Shear Ratio 
0.82 F 
0.80 F 
vl 

0.78 F 
0.76 F 
0.747 
0-79.50 0.25 0.30 0.35 0.40 0.45 


FIG. 25. A test of the impact of alternative redshift analysis choices 
on the inferred cosmology from 3x2pt. We compare the fiducial 
3x2pt analysis (black) to an analysis where we marginalize over the 
ensemble of n(z) realizations directly via Hyperrank instead of their 
effective mean redshifts (blue) and to an analysis where we remove 
the shear-ratio data (orange). 


Fid. 3x2pt 
0.84r ADM Nonlinear Bias 
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FIG. 26. A test of the impact of alternative analysis choices on 


the inferred cosmology from 3 x2pt. We compare the fiducial 3x2pt 
analysis (black) to an analysis where we marginalize over a nonlinear 
bias model using smaller scales in y and w(0) in blue, an analysis 
that marginalizes over free lens magnification bias parameters in or- 
ange, and an analysis that uses the NLA IA model in green. 
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TABLE V. The MagLim galaxy bias constraints (mean with 68% CL) from the fiducial linear bias analysis and the nonlinear bias analysis. 
Due to a nontrivial second peak in the posterior of the b{, we also show in parentheses the ID marginalized peak value. 


by m b3 ba b? b2 bZ b; 
Fid. Linear Bias 1.49019 1.697011 1.917073 1.79701] — - - - - 
Nonlinear Bias 1.447702 1.604810 1.854510 1.74% air O.12* 1g 0.067026 —U.06* 627 0.401982 (0.71* 0255 
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FIG. 27. The validation of the model and inference pipeline on 


the Buzzard simulation suite, where the true cosmology is indicated 
by the cross. An analysis of a synthetic data vector (black) at the 
true Buzzard cosmology based on the true redshift distributions with 
fixed shear and photo-z parameters is compared to a full analysis of 
the mean data vector of 18 simulation realizations (blue) including 
all nuisance parameters and n(z) distributions inferred in the same 
way we do using the real survey data. 


Appendix E: Alternative analysis choices and robustness 


In addition to the validation described in the main text of 
this work and the associated papers described in App. |A| we 
discuss several other analysis modifications in this Appendix 
that test the robustness of our result. 


1. Photometric redshifts 


In Fig. [25] we show the fiducial 3x 2pt analysis compared 
to two variations in how we use information related to the 
photo-zs of our source sample. First, we simply remove the 
shear-ratio part of the data vector, which helps to constrain the 
photo-z parameters, but does not contribute directly to con- 
straining cosmological parameters. Second, we marginalize 
over the full ensemble of n(z) realizations. Before unblind- 
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FIG. 28. A comparison of the fiducial 3 x 2pt analysis (black solid) 
with one that fixes the neutrino mass density to fixed is minimum 
value (black dashed). We make a similar comparison for the Planck 
CMB data (green solid and dashed). 


ing, we decided to simplify the way we use the redshift infor- 
mation for the sources and marginalize over the more typical 
set of parameters that encode shifts in the mean redshift of the 
average n(z) from the ensembles. This choice was driven by 
the additional computational expense of marginalizing over 
the n(z) realizations directly, but was demonstrated in sim- 
ulated analyses to produce consistent cosmological results at 
DES Y3 precision. To test this, we run a single chain to show 
the potential impact of this choice on the real data. While 
all three results are consistent, with a shift of 0.530 between 
Hyperrank and marginalizing over Az, it is likely that future 
analyses will need to marginalize over the full ensemble of 
n(z) realizations directly using something like the Hyperrank 
process. Finally, we utilize the SOMPZ framework to red- 
erive the redshift distributions of the magnitude limited sam- 
ple [126], instead of relying on the fiducial DNF redshifts. 
The 3 x 2pt cosmology inferred from this alternate redshift 
distribution is fully consistent with the fiducial result. 
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TABLE VI. Summary of internal goodness-of-fit tests with the posterior predictive distribution for analyses with the redMaGiC and MagLim 
samples. The first and second columns indicate the subset of data d considered. Realizations of d are generated for model parameters drawn 
from its own posterior and compared to actual observations of d, following the method developed in Ref. [132]. The third column shows the 
(calibrated) probability-to-exceed for the full data test used for comparison. The other columns indicate the probability-to-exceed for different 


subsets used for comparison. 


Goodness-of-fit test d 


p(d|d) p(&x|d) p(w + dd) p(é+ld) p(§-ld) pCvild) p(wld) 


MagLim 

3x 2pt €+,7%,w 0.023 0.209 0.209 0.243 0.353 0.009 0.296 
2x 2pt ^«,w 0.019 E 0.019 - - 0.013 0.356 
Cosmic shear E+ 0.218 0.218 - 0.220 | 0.359 - - 
redMaGiC 

3x2pt E+, yt, w 0.005 0.242 0.001 0.260 0.410 0.029 0.041 
2x2pt ^«,w 0.038 - 0.038 - - 0.187 0.072 
Cosmic shear E+ 0.255 0.255 - 0.253 0.397 B Š 


TABLE VII. Summary of internal consistency tests with the posterior 
predictive distribution for analyses with the redMaGiC and MagLim 
samples. We consider the consistency between different two-point 
functions, as listed in the first column. The probability-to-exceed 
p(a|b) is obtained by generating realizations of a for model parame- 
ters drawn from the posterior of b, and comparing these realizations 
to the actual observations of a. The third and fourth columns indicate 
the probability-to-exceed for the redMaGiC and MagLim samples re- 
spectively. 


Consistency test p redMaGiC MagLim 
Cosmic shear vs galaxy-galaxy p(€+|%,w) 0.019 0.297 
lensing and clustering 

galaxy-galaxy lensing vs cos- p(y%|€4,w) 0.016 0.004 
mic shear and clustering 

Clustering vs cosmic shear and p(w|&+,7%) 0.008 0.526 


galaxy-galaxy lensing 


2. Nonlinear bias modelling 


The fiducial analysis of galaxy clustering and galaxy— 
galaxy lensing assumes a linear galaxy bias model, which 
requires removing significant small-scale information in our 
data vector that exhibits substantial nonlinear behavior. This 
keeps the analysis simpler with fewer free parameters, but po- 
tentially wastes this additional information to further constrain 
cosmology on small scales. In our fiducial model, we only 
utilize large scales above 8 ^ ! Mpc for w(0) and 6 h^ Mpc 
for y+ measurements. To remedy this, we have developed a 
nonlinear galaxy bias model and analysis that can utilize this 
smaller scale information from the lens sample at the expense 
of marginalizing over additional bias parameters. This is sim- 
ilar to going to smaller scales in cosmic shear and having to 
marginalize over additional baryonic effect freedom. In both 
cases, the improvements are limited by needing to simultane- 
ously constrain these additional parameters. 

Our model is a hybrid 1-Loop effective field theory model, 
having five free parameters, as detailed in Refs. 
[141]. In Ref. we validated this model using 3D 
correlation function measurements in redMaGiC and MagLim 
mock catalogs and find that this nonlinear bias model agrees 
with significantly higher signal-to-noise measurements at bet- 


ter than 2% above scales of 4 h~'Mpc. For cosmological in- 
ferences we fix three of these parameters based on theoretical 
considerations and validation on these mock catalogs in or- 
der to reduce parameter degeneracies and potential parameter 
volume effects. We have validated the analysis using multiple 
buzzard simulation realizations([119][140] [141] (see Sec. 
for details), finding that this model gives a cosmological bias 
of less than 0.30 in the Q,,—53 plane for both the w(0) + y: 
and 3x2pt probes. 

In Fig. 26] we show the cosmological constraints when ap- 
plying this model to our data. The galaxy bias constraints 
are summarized in Table[V] We find cosmological constraints 
in wCDM to be consistent with the fiducial analysis. We 
find a similar improvement in constraining power for ACDM 
with this nonlinear bias analysis as with the ACDM-optimized 
analysis that instead incorporates additional small scale infor- 
mation in cosmic shear. With DES Year 6 data, it is possible 
we could better optimize these choices of how to utilize small- 
scale information in ACDM analyses to substantially improve 
constraining power. 


3. Lens magnification bias parameters 


In the fiducial analysis, we fix the lens magnification co- 
efficients C7 to values derived in Ref. [139]. In Fig. [26] we 
show the impact of freeing these coefficients with a wide flat 
prior between —6 and 10 in each lens bin. We find the cosmo- 
logical constraints in ACDM to be consistent with the fiducial 
analysis. Further investigations into the lens magnification co- 
efficients and their impact on cosmological constraints can be 


found in Ref. [139]. 


4. Intrinsic alignment models 


In the fiducial analysis, we use the full TATT intrinsic align- 
ment model. In Fig. [26] we show the impact of limiting the 
intrinsic alignment model to the more commonly used NLA 
model with free redshift power-law evolution. This is the 
same as fixing the parameters A5 = bra = 0 in the TATT 
model. We find consistent cosmological parameters, with a 
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TABLE VIII. A comparison of metrics testing the consistency of independent data sets within the ACDM model. We show results in units of 
c, with a derived probability-to-exceed p shown in parentheses. Overall, we find no significant (defined as p < 0.01) evidence of disagreement 
between the DES 3x 2pt, external low-redshift, or Planck CMB data sets. The details of the metrics and derivation of p are described fully in 


Ref. [133]. 

Data set 1 Data set 2 Parameter difference Suspiciousness Eigentension Qupw  ln(Evidence ratio) 
DES 3x2pt Planck CMB 1.50 (0.13) (0.7 + 0.1)o (0.48 + 0.08) 1.20 (0.23) 1.60 (0.10) 5.90.5 
DES 3x2pt Ext. Low-z 1.30 (0.11) (1.6 + 0.3)e (0.13 + 0.07) 0.90 (0.40) 1.40 (0.11) 1.6 + 0.4 
Ext. Low-z Planck CMB 1.1c (0.27) - 1.50 (0.13) 0.66 (0.54) - 

DES 3x2pt + Ext. Low-z Planck CMB 0.90 (0.34) - 1.20 (0.25) 0.8c (0.42) - 

DES ACDM-Opt. 3x2pt Planck CMB 1.9c (0.06) - 1.30 (0.21) - - 

DES Hyperrank 3x2pt Planck CMB 2.20 (0.02) - 2.00 (0.04) - - 

DES redMaGiC 3x2pt Planck CMB 1.50 (0.13) - 1.00 (0.30) - - 


gain of only 13% relative to TATT in the og—Q,,, parameter 
plane. The use of NLA was not motivated a priori as the fidu- 
cial IA model, but due to the IA amplitude parameters being 
constrained to be small in TATT, Ag = bra = 0 is not a poor 
approximation. 


5. Simulation validation tests 


We reproduce the fiducial analysis on a suite of 18 Buzzard 
simulations described in Sec.[IV DJ This is shown in Fig. 27] 
where the true cosmology is indicated by the cross. We com- 
pare two simulated analysis. The first analysis uses a syn- 
thetic, noiseless data vector based on the true n(z) from Buz- 
zard and without marginalizing over shear or photo-z parame- 
ters. The second reproduces the full analysis on the mean data 
vector of all 18 simulation realizations, which marginalizes 
over all nuisance parameters and uses an photo-z n(z) that is 
inferred from the same process we apply to the real survey 
data. We find that the two simulated analyses agree very well 
with each other and with the true cosmology for each simu- 
lated data vector. 


6. Neutrino mass 


Finally, we also show a version of the analysis that fixes 
the neutrino mass density at the minimum allowed mass in 
Fig.28] comparing to a similar fixed neutrino mass density of 
the Planck CMB data. 


Appendix F: Details of Internal and External Consistency 


Results of the final internal consistency tests are reported in 
Tables [VT] and Table [VI|reports results of the goodness- 
of-fit tests where realizations of a subset of the data vector are 
generated for model parameters drawn from its own posterior. 


We show results for analyses using the MagLim lens sample 
and for redMaGiC. Table [VII] reports extra consistency tests 
between the two-point functions that were not included as un- 
blinding criteria. In these consistency tests, the data vector 
is split in two, and one part is used to predict the other. For 
instance, the first row of Table [VII] shows the result of a test 
where realizations of cosmic shear are generated for model 
parameters drawn from the 2 x 2pt posterior, and compared 
to observations of cosmic shear. All pre-defined unblinding 
requirements were met for these internal data combinations, 
and both 3x2pt and the cosmologically-constraining subsets 
of the data for both lens samples show consistency with each 
other and with ACDM. However, there persists evidence of 
potential issues with the consistency of some parts of the com- 
bination of galaxy clustering and galaxy-galaxy lensing rela- 
tive to the best-fit model of the full 3x2pt or cosmic shear 
data. In particular, for redMaGiC this seems slightly more 
likely to be sourced primarily from the clustering data. 

We report detailed consistency metrics between external 
data set pairs in Table The details of the metrics and 
derivation of probability-to-exceed p are described fully in 
Ref. (133]. Overall, we find no significant (defined as p < 
0.01) evidence of disagreement between the DES 3x2pt, 
external low-redshift, or Planck CMB data sets shown in 
Figs. [14] & [15] We show both a method measuring parameter 
differences and a method based on the evidence (‘Suspicious- 
ness’) for the two primary comparisons of DES 3x 2pt with 
external data sets. The evidence ratio is also shown, and all 
three metrics give qualitatively consistent results. 

All data set combinations are largely in agreement, accord- 
ing to all tension metrics. Some differences between methods 
are likely due to non-Gaussianity of both posteriors involved, 
as well as the different approaches that the methods employ in 
treating the impact of priors on tension quantifications, but in 
all cases are below a fraction of a sigma. We also note that the 
Eigentension result for the case of DES 3 x 2pt vs complemen- 
tary external low-redshift probes may not be capturing the full 
tension; the two data sets have comparable constraining pow- 
ers, for which case Eigentension is not optimized. 


